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We consider theoretically the electron-electron interaction induced exchange- 
correlation effects in the lowest subband of a quasi-one-dimensional GaAs quantum 



J> . wire structures. We calculate, within the leading order dynamical screening approxi- 

mation (i.e. the so-called GW approximation of the electron gas theory), the electron 
self-energy, spectral function, momentum distribution function, inelastic scattering 

m ; 

ON . rate, band gap renormalization, and, the many-body renormalization factor both at 



zero and finite temperatures, and, both with and without impurity scattering effects. 
We also calculate the effects of finite temperatures and finite impurity scattering on 



C . 

O ■ the many-body properties. We propose the possibility of a hot-electron transistor 

O 

device with a large negative differential resistance which is based on the sudden on- 
set of plasmon emission by energetic ballistic electrons in one dimension. The issue 
of the existence or non-existence of the Fermi surface among the interacting one- 
dimensional quantum wire electrons is critically discussed based on our numerical 
results. 
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I. INTRODUCTION 



Recently, there has been a great deal of interest |IJ in ultranarrow confined semiconductor 
systems, called quantum wire structures, where the motion of the electrons is essentially 
restricted to be one- dimensional. These systems are usually fabricated on high-quality two- 
dimensional electron systems by confining the electrons in one of the remaining free directions 
through ultrafine nanolithographic patterning. ^IM|3|| In addition, other methods for 
fabricating wires, such as growth of wires on tilted superlattices and vicinal substrates 0] 
have also been reported. These developments have opened up possibilities for novel and 
exciting experiments, and stimulated considerable theoretical activity. [§,^|, |lO|)|TTi|l2i|l3| At 
present, quantum wires with active widths as small as approximately 300 A and of negligible 
thickness have been fabricated, with further reductions in size and improvement in quality in 
the offing. Quantum wires have generated much interest both for the potential for practical 
applications in solid-state devices such as high-speed transistors, efficient photodetectors 
and lasers, fll^ , |T5|| and because they have afforded us for the first time an experimental 
opportunity to study real one-dimensional Fermi gases in a relatively controlled manner 
(much in the same way that, in the past two decades or so, semiconductor inversion layers, 
hetero junctions and quantum wells have generated considerable activity in pure and applied 
research of the two-dimensional electron gas). Thus, both from the fundamental and applied 
physics viewpoints, there is interest in understanding the electronic properties of quasi-one- 
dimensional quantum wires. 

Much of the fundamental theoretical understanding of electrons in one- dimensional sys- 



tems have come from work on the Tomonaga-Luttinger model. ||i£^|17| , |18f The Tomonaga- 
Luttinger model makes some drastic simplifying assumptions which allow one to solve the 
interacting problem completely. One surprising result that is obtained from the solution 
of this model is that even the smallest interaction results in a disappearance of the Fermi 
surface, leading to a system which is non-Fermi liquid (in the sense that the elementary ex- 
citations are very different from those of the non-interacting system). Therefore, one would 
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expect that the experimental properties of the semiconductor quantum wires should be very 
different from any predictions based on the assumption that the one-dimensional electron 
gases are Fermi liquids. 

Recently, quantum wires have been fabricated |4]|^,fJ in which the one-dimensional quan- 
tum limit has been attained, in the sense that only one quantum subband is populated by 
the electrons in the quantum wire, so that a one- dimensional interacting Fermi gas model is 
valid. Contrary to what one might expect based on the well-established theoretical results, 
the experimental data obtained from these experiments on these structures, such as Raman 
scattering and photoluminescence, can be successfully explained on the basis of standard 
Fermi liquid theory. This fact is also mirrored by the success of other theories, such as the 
Landauer-Biittiker formula, [|19j which treat electrons in mesoscopic quasi-one-dimensional 
systems as noninteracting entities. Thus we are posed with the question, "Why are experi- 
mental results from real quantum wires explicable on the basis of Fermi liquid theories?" This 
is one of the points we attempt to shed light on in this paper. We show that the interplay 
between plasmons and impurities play an important role in determining the physical char- 
acteristics of the system. In particular, we show that the physical mechanism responsible 
for the disappearance of a Fermi surface in a clean system is the virtual plasmon emissions, 
and that in an impure system these virtual plasmons emissions are damped, leading to a 
restoration of the Fermi surface. (The plasmons in our model are equivalent to the Tomonaga 
bosons of the Tomonaga-Luttinger model.) 

We also present calculations for experimentally measurable properties for quantum wires. 
For the quantum wire to be in the quantum limit, the doping must necessarily be low 
< 10 6 cm _1 (otherwise, higher subbands are filled) and hence the Fermi energy and other 
energy scales pertinent to the system are also small. Therefore, temperature effects are 
important, even at relatively low temperatures. We have therefore presented the calculations 
for both T = and finite temperatures. To perform the T ^ calculations, we have 
developed a formalism for the efficient calculation of the finite-temperature self-energy (valid 
for arbitrary dimensionality), which we describe in this paper. 
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The many body properties of one-dimensional systems have unique properties, one of 
which is the singular behavior of the mean free path as a function of electron energy, which 
occurs when electrons cross the plasmon emission threshold. We discuss the possibility that 
this singular behavior could in principle be applied to produce a device with large negative 
differential resistivity. Finally, since the the one- dimensional plasmons are prominent pro- 
tagonists in the physics of interacting one-dimensional system, we conclude the paper with 
a discussion of the plasmons in the presence of impurities and at finite temperature. 

This paper is organized as follows. In Sec. |T|, we study the many-body properties of 



quantum wires at T = 0. In Sec. |T|, we develop the formalism (for arbitrary dimension- 
ality) necessary for efficient calculation of the self-energy at finite temperatures, which we 
apply to calculate the many-body properties of a quantum wire. In Sec. |TV|, we discuss the 
possible device applications of quantum wires. In Sec. [V|, we discuss the plasmon spectrum 
of quantum wires, in clean and dirty systems, and at finite temperature. Sec. [V| contains our 
summary and conclusions. Some of the results in this papers have been published previously 
in brief communications. HJT^JTH] Also note that in this paper, h = 1 and k& = 1. 



II. MANY-BODY PROPERTIES OF A QUANTUM WIRE AT T = 

In this section, we present calculations of the many-body properties of a one-dimensional 
quantum wire, using the so-called GW approximation. We begin with a brief review of the 
formalism, and then we describe how we put effects of impurities into the system. We then 
give results and conclude this section with a discussion. 

We assume that the quantum wire is formed by confining the electrons to the two- 
dimensional plane x—y plane (by, say, modulation doping an AlGaAs-GaAs-AlGaAs het- 
erostructure), and then further confining the electrons along the ^/-direction. We assume 
that the confinement in the z-direction is much stronger than the confinement in the y- 
direction, which is a reasonable approximation because, at this stage, the technology for 
confining electrons to two-dimensions is much more advanced the technology for confining 



4 



the electrons along an additional direction. For example, by modulation doping, the elec- 
trons can be confined in the z-direction on the order of less than 100 A, whereas currently 
the best confinement in the y-direction is approximately 300 A, leading to at least an order 
of magnitude difference in the energy level spacings of the y- and z- directions. Thus, in 
this paper, we assume that the electron gas has zero thickness in the z-direction, but has a 
finite width in the y-direction. (Note that it is necessary to assume a nonzero width so that 
the one-dimensional Coulomb interaction is finite. |]20[]) The confinement of the electrons in 
the y-z plane leads to the quantization of energy levels into different subbands, depending 
on what the wavefunction in the y-z plane is. We assume that the energy separation be- 
tween the lowest energy and higher energy subbands is so much larger than all other energy 
scales in the problem that the higher subbands can be ignored. Therefore, all quantities in 
the main part of this paper refer to the lowest energy subband. In Appendix A, we give 
the formalism in the case of multiple subbands and briefly discuss why, for processes with 
energies much less than the subband energy separation, the higher subbands have negligible 
contribution. 



A. The Hamiltonian 

The Hamiltonian of the system is given by H = H + H c - C where Ho is the bare kinetic 
energy term and H e _ e is the electron-electron interaction term, given by 

ka 

(2.1) 

H e-e = ^ J2 12 V M) d k+q,a A* -q,a> CfcV C k a, 
kk'q era' 

Here c ka (c^) are fermion creation (annihilation) operators for for states with wavevector 
k in the x-direction and electron spin a, £(k) is the is the kinetic energy relative to the 
chemical potential fi, and L is the length of the system. The V c (q) are the matrix elements 
of the Coulomb interaction (screened by the static lattice dielectric constant eo), 
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where <fi(y) is the transverse wavefunction and 



2e 2 



v(q,y-y') = — K Q (\q(y - y')\) 



(2.3) 



is the one-dimensional Fourier transformation of the Coulomb interaction. [^,^T|] Here, 



Kn{x) is the zeroth-order modified Bessel function of the second kind. 22 



For an external confining potential that is parabolic, the self-consistent confining po- 
tential (that is, the external confining potential plus the self-consistent Hartree term) gives 
approximately a square well potential. [23] This turns out to be a reasonable approxima- 
tion for modeling quantum wire confinement. Therefore, we approximate the self-consistent 
potential by a square well with infinite barriers at y = —a/2 to y = a/2, so that 

r 2 



cos 



(^), if -a/2 < y < a/2; 



<f)(y) = < 



(2.4) 



0, otherwise; 
and the Coulomb matrix element in this square-well case is 



2e 2 

VM) = — 



2 fa/2 



a/2 



dy I dy' K Q (\q(y - y')\) cos 2 (— ) cos 2 (— 

a/2 J-a/2 K a 7 K a 



(2.5) 



2e 2 f 1 r 3 

— / dx K (\qa\x) (2 — (1 — x) cos(27ra;) + — sin(27ra) 
eo Jo 1 2ti 



whose asymptotic forms are 

f 37re 2 



VM) 



e \qa\ 



for \qa\ — > oo; 



(2.6) 



2e 2 

[K {\qa\) + 1.9726917...], for \qa\ -> 0. 



6 



Note that Kq(x) ~ — ln(x) as x — > 0, so that the Coulomb matrix element diverges as 
— 2e 2 ln(|ga|)/eo small q. To speed up numerical computation, we fitted V c (q) with a function 
which extrapolated between the q — > and q —>■ oo forms and which deviated from the true 
function by no more than 0.6%. 

Throughout this paper, we assume that the band for the electrons in the unconfined 
^-direction is parabolic, so that = h 2 k 2 /(2m*) — ji, where m* is the effective electron 
band mass. This is an excellent assumption for GaAs quantum wires because the relevant 
electron densities are so low (kp ~ 10 6 cm -1 ) that the nonparabolicity in the band for k ^ fcp 
is insignificant. 

B. The electron self-energy H(k,oj) within the GW and RPA approximation 

A quantity which provides a substantial amount of information on an interacting many- 
electron system is the electron Green's function G(k,u), or equivalently, the electron self- 
energy, H(k,u) = G' Q 1 (k,uj) — G~ l (k,uj) (where Go(k,u) is the bare noninteracting Green's 
function). The self-energy is roughly the correction to the noninteracting electron single- 
particle energy due to the interactions. Once G(k,u) or Ti(k,u) is known, the one-electron 
properties of a system, such as the spectral density function, electron distribution function 
and band-gap renormalizations, can be calculated. 

Generally, it is impossible to calculate the u) exactly for interacting systems, and one 
must resort to various approximation schemes. However, in one dimension, the self-energy of 
the Tomonaga-Luttinger model, which uses somewhat artificial assumptions (specifically (1) 
two completely linearly dispersing bands of electrons with an infinite bandwidth, populated 
by an infinite density of electrons and (2) a short-ranged interaction), can be calculated 
exactly. [fTl| , |I7|1 As alluded to in the introduction, this models predicts that any interaction, 
no matter how small, will drive the system away from Fermi liquid behavior. However, the 
connection between the exactly solvable but artificial Tomonaga-Luttinger model and real 
semiconductor quantum wires is somewhat tenuous, and there is a need for calculations on 
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more realistic models which can be compared to experiment. The model we employ assumes 
that the electron dispersion in the free direction of the quantum wire is parabolic, and 
the electrons interact via the exact Coulomb interaction. Unfortunately, with these more 
realistic assumptions, the problem is no longer exactly soluble and therefore, as in higher 
dimensions, u) can only be calculated approximately. In this paper, we use the so-called 
GW approximation to calculate the self-energy, |24| where the Green's function is dressed 



by the dynamically screened Coulomb interaction, which is calculated within the random 
phase approximation (RPA). The Feynman diagram for the self-energy and the screened 
interaction in these approximations are shown in Fig. |l](a) and [l|(b). This approximation, 
which sums up the largest diagrams in the self-energy perturbation series, has long been 
employed to successfully calculate properties of three and two dimensional electron systems. 
|25| , |2~6|j We ignore diagrams which are higher order in the screened interaction, such the 
diagram shown in Fig. 0(c), because we expect them to have a smaller contribution. In fact, 
in Appendix B we show that at k = kp and uj = and for a short-ranged interaction, 
Fig. [TJ(c) gives zero contribution. 

The self-energy within the GW approximation at T = is 

E(k, u) = if^r^ W(q, uj')G (k -q,u-cS). (2.7) 

where, G (k, uj) — [w — £/. +iO + sgn(/c — kp)]' 1 and W(q, uj) is the screened Coulomb inter- 
action, which given by 

W(q,uj) = ^- (2.8) 

Here, e(q, uj) is the dielectric function, which describes the screening properties of the elec- 
tron gas. The self-energy can be separated into the frequency-independent exchange and 
correlation parts 

E(k,cu) = Z cx (k) + Z COT (k,u). (2.9) 
The exchange part is given by Jl8| 



— n F (k + q)V c (q), 



(2.10) 



where np{k + q) — 6{k, F — \k + q\) is the Fermi function at T = 0, and £ cor (fc,u;) is defined 
to be the part of not included in S ex (/c). In the GW approximation, the H COT (k,uj) 



can be written in the line and pole decomposition [18.27 



S cor (/c, uj) = Si ine (A;, u) + E P oie(A;, uj) 

„ , r°° dq„ , , r°° du' 

Ei ine (/C, uj) = - 



00 Z7T J -00 Z7T [qk+q — UJ) 



1 



J pole 



(k,u) 



00 |- 
-00 2tt 



(u-S k+q )-0(-S k+q )]V c (q)( 



e(q,iu') 
1 



- 1 



(2.11a) 
(2.11b) 

(2.11c) 



The Ei ine (/c, a;) is completely real because e(q,iuo') is real and even with respect to to'. Thus, 
ImfSpoie] gives the total contribution to the imaginary part of the self-energy. In evaluating 
the self-energies, we assume that the frequency u has an infinitesimal positive imaginary 
part; this gives us the retarded self-energy. |Tj| (Unless otherwise stated, all real frequency 
self-energies in this paper are retarded.) 

Within the RPA, e(q,u) is given by e(q, uj) = 1 — V c (g)rio(g, u>), where Uo^,^) is the 
irreducible polarizability, given diagrammatically by the particle-hole bubble. In a pure 
system, PJ20IJ21 



n (g,z) 



rn 



In 



nq 



z 2 — ((q 2 /2m) — qv F y 



(2.12) 



-z 2 — ((q 2 /2m) + qv F ) 2 - 
where the principal value of logarithm (|Im[ln] | < it) should be taken. In evaluating Tlo(q, uj) 
for real frequencies, the limit z = uj + i0 + should be taken. The RPA form of e(q, uj) has 



recently been shown ]29| to exactly reproduce the plasmon dispersion of one-dimensional 
systems, in the sense that the long wavelength RPA correctly gives the exact Tomonaga 
boson dispersion of the Tomonaga-Luttinger model. 

The Ho(q, uj) given in Eq. ( |2.12j ) assumes that the system is free from any defects. While 
modulation doping techniques result in quantum wires are generally of high quality, impu- 
rities and other imperfections always exist. In general, the effect of these impurities is to 
cause the electrons to diffuse instead of moving ballistically. This diffusive motion of the 



electrons affects the polarizability (since it is the response of the electrons to an applied 
potential) and hence electronic screening, generally tending to make screening less effective 
since electrons cannot move as quickly to their screening positions. Impurity effects are 
usually introduced diagrammatically into the screening in the RPA by including impurity 
ladder diagrams into the electron-hole bubble, [[HJ which yields a polarizability H y (q,uj) 
which has a diffusive regime for small q and uo. Since the exact expression for n 7 (g, u) 
within this diagramatic approach is complicated, we use a particle-conserving expression, 



given by Mermin, |3T| which captures the essential physics of impurity collisions on II 7 (g, u). 
For an impurity scattering rate 7, 

U ^'^-o; + i 7 [no(g,a; + z 7 )/n ( g) 0)]- 
The form given in Eq. (|2.13|) has the correct diffusive behavior at small q and lu, which is 

absent in the n (g, lu) for the pure case. In particular, for q,u — ► 0, 

rr / x 2m e Dq 2 

irk F n Dq 1 + iuj 

(2.14) 

9 

U (q,u)^^—r (for q < w/uf), 

where D = is the one-dimensional diffusion constant. The modification of the behav- 
ior of the polarizability due to impurities has important consequences for long-wavelength 
plasmons, which in turn affect the many-body properties of the one-dimensional electron 
systems significantly. 

For frequencies on the imaginary axis (as is the case in the integrand of the line component 
of the self-energy) , the polarizability is given by 

o / ■ i\ (|w'| + 7)no(g,i|u/|+i7) , . 

n 7 (g, iu ) = — — — — - — , (2.15) 



C. Calculation and Results 



Using the expressions given above, we calculate the self-energy for a quantum wire for 
both pure and impure cases. We briefly discuss some subtleties involved in the numerical 
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calculation, and then we show our results. 

Both the S ex (fc) and Ei ine (/c, to) are comparatively straightforward to evaluate numeri- 
cally. However, evaluation of E po i e (fc, to) for the clean case (7 = 0) is complicated by the 
singularities present in the integrand. 

For Im[E po ie(fc, a;)], in addition to the contribution from single-particle excitations 
(Im[e] 7^ 0), there are also contributions from the plasmon excitations (Re[e] = and 
|Im[e] I = + ), which produce ^-functions in the integrand of Im[S po i e (/c, uj)] (see Fig. §). Care 
must be taken to pick up this ^-function in a numerical integration. To numerically evaluate 
the plasmon contribution to Im[£], we performed a search along q for Re[e(g, £k+ q — &)] = 0, 
in the regions Im[n (g, (,k+q — <^)] = 0. When a zero of Re[e] was found, we evaluated the 
integral by two methods: (1) introducing a small damping term 7 in e and numerically 
integrating around the Re[e] = region and (2) numerically evaluating de(q,£,k+ q — u>)/dq, 
which is proportional to the weight of the 5-function. Both methods gave the same results. 
Since the danger of missing the plasmon peak in the case of small 7 is always present, we 
routinely searched for the Re[e] = even when 7 7^ 0. 

In the integrand for Re[£ po i e (&, uj)], at the q where Re[e] = 0, there is a (g — go) 1 
principal part divergence, which also can cause numerical problems if it is not handled 
carefully. We treat this the problem either by cutting off the integrand around go when its 
absolute value exceeded a certain value or by putting in a small damping term 7, and paying 
special attention to the integration around go- Again, both methods gave the same results. 



From £(£;, oj), we can obtain the single-particle spectral function [0,^2 



m u) = 2|Im[£(^)| (2 16) 

l ' ] (cu-e fc -Re[E(A;,u;)]) 2 + (Im[S(fc,^]) 2 . K ' ] 

(Note that Im[E] is assumed to contain an infinitesimal negative part.) The spectral func- 
tion A(k, uj) can roughly be interpreted as the probability density of the different energy 
eigenstates required to make up a particular /c-state. It satisfies the sum rule 

±-A{k,u>) = l, (2.17) 
00 Ztt 
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which is generally satisfied to within less than a percent in all our numerical calculations. 
We pay particular attention to the self-energy and spectral function at k — kp, since the 
behavior of of these functions at the Fermi surface determine the low energy properties of 
the system. The exact solution of the Luttinger model indicates that interactions produce a 
non-Fermi liquid, the so-called Luttinger liquid. We reproduce this result in clean systems, 
and show that virtual plasmon emission is responsible for this behavior. We then show that 
impurities suppress virtual plasmon emission, and the Fermi surface reappears. 

A system is a Fermi liquid if it possesses a Fermi surface (i.e., a discontinuity in n&) 
whose presence is indicated by a 5-function in A(kp, uj) at uj = 0. The existence of a S(uj) in 
A(kp,u) depends crucially on the behavior of lm[T,(k F ,u)] as uj — > 0. If |Im[£(fcp,u;)]| goes 
to zero faster than M, and 



u - 6=f - Re[E(k F ,u)} m Z f 1 uj, 



where Zp is a constant called the renormalization factor, |u| then as uj 
|Im[E]|, and hence 

2e 

A{k F , uj) lim 



(2.18) 



0, \Re\L{k F ,uj)\ > 



*-o Z~ f x uj 2 + e 2 
2% Zp 8(u>). 



(2.19) 



Thus, if | Im uj)]\ goes to zero faster than |a;|, then there is a discontinuity of magnitude 
Zp in rifc at kp. Since Zp is the linear coefficient of the expansion in uj of the left hand side 
of Eq. ( |2.18|) , it is given by 

-dRe[E(k F ,u)]- 



duo 



u=0 



^2.20) 



In contrast, if |Im[2D(A; j p, cu)] | goes to zero slower than then the spectral function 
A(kp,u) does not take the form Eq. ( |2.19| ) as uj —>■ 0, and hence there is no 5-function in 



A(kp,u), implying that the system is not a Fermi liquid. ||l6| , |lT| , |33f In general, interacting 
three-dimensional systems with and without disorder, and, interacting pure two-dimensional 



systems are Fermi liquids. |3^ , p5| , |36|] Through a study of Im[£(/jp, uj)], we show that in one 
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dimension within the GW approximation, the system is not a Fermi liquid in the absence 
of impurity scattering, but becomes a Fermi liquid in the presence of scattering. 

The imaginary part of the self-energy in the GW approximation, from Eq. (p.llcQ , can 
be written as 



1 



OC 



lm[E(k, uj)] = -j_^ [9{u - £ fc+g ) - e(-Zk +q )\ P(q, £*+, ~ u), (2.21) 

where P(q,E) = 2V c (q)lm[e~ 1 (q, E)] which is the Born approximation transition rate for a 
momentum change q and energy change E. [^7] Thus, Im[E(/c, u)] can be interpreted as the 



total scattering rate for an electron with momentum k and energy uj. By setting u = the 
actual electron kinetic energy, and one obtains the total Born approximation scattering rate. 
For uj ^ £k, Jm\E(k,u)] corresponds to the total scattering rate due to virtual transitions, 
since the transitions violate energy conservation. 

At low energies in two and three dimensions, single-particle scattering is far more im- 
portant than plasmon scattering because the single-particle excitation spectrum is gapless 
and the phase-space available for single-particle scattering extends around the entire Fermi 
surface, whereas the plasmon dispersion either rises quickly or has a gap at q = 0. Therefore, 
for small uj, the major contribution to Im[£(&F, uj)] in two and three dimensions comes from 
virtual single particle excitations. In contrast, in one dimension, the single-particle excitation 
spectrum has a gap except at \q\ = 0,2kp, and the phase-space available for single-particle 
scattering is severely restricted, while the plasmon dispersion is gapless at q = 0. Hence, in 
one dimension, lm[Tj(kp, uj)] at small uj is dominated not by virtual single-particle excitations 
but by the virtual excitation of plasmons. It is the domination of the plasmon excitations 
at low energies that give interacting one- dimensional systems their unique non-Fermi liquid 
behavior. 

In the case of a clean quantum wire (7 = 0), as uj — > 0, the contribution to Im[£(fcp, uj)] of 
the low energy plasmon excitations goes as |a>| | ln( |o; | ) | 1//2 , whereas the contribution of the 
single-particle excitations goes as u;/(ln M) 2 . Thus the plasmon contribution dominates, 
and since Im[E(&F, u)] goes to zero slower than |u;|, the Fermi surface does not exist (in 
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agreement with Luttinger liquid theory). The Fermi surface is smeared out to the extent 
that a momentum-space discontinuity in the distribution function does not exist because 
of the ease with which the particles at the Fermi surface can emit virtual plasmons. In 
contrast, in two dimensions, the contribution to Im[E(A;.F, u;)] of low energy plasmons goes 
as uj 2 , whereas that of the single-particle excitations goes as u 2 \ ln(|u;|)| and therefore in 
the two-dimensional case, the single-particle excitations dominate, and since the self-energy 
goes to zero faster than \u\, the two-dimensional interacting electron gas is a Fermi liquid. 

The inclusion of impurity scattering causes the electrons to diffuse at long wavelengths 
which damps out the plasmons at small q (see section 0). Hence, the plasmon contribution to 
Im[S(A;jr, (jS)\ at small \oj\ is removed, and the entire contribution comes from single-particle 
excitations. In one dimension, it modifies the behavior to |Im[£(/cp, ~ uo 2 \ ln(|u;|)| 3 as 
\uj\ — > 0, which implies that the Fermi surface is restored. This result indicates that the 
Fermi surface is resurrected in dirty systems because the low energy virtual plasmon emission 
responsible for its destruction in clean systems has been suppressed by impurity scattering. 
In two dimensions, however, the situation is reversed. The plasmon contribution was not 
dominant in the first place, so its removal is not significant. However, inclusion of impurity 
scattering enhances [^SJ the single particle scattering rate (because of phase space effects) 
and it results in a |Im[£(A;_F,a;)]| ~ This agrees with the notion that an interacting dis- 
ordered two-dimensional system is not a Fermi liquid. |38| Table | summarizes the behavior 



of [£(fcp,u;)] resulting from various contributions in both the clean and dirty cases, in both 
one and two dimensions. The details of the calculations are given in Appendices C and D. 

In Fig. ^ we show the real and the imaginary parts of the self-energy for k = (bandedge) 
and k = kp (Fermi surface), for E = 0,Ep, for kpa = 0.9 and r s = 4e 2 /(jreoVp) = 1.4. 
f39|j In the k = kp, clean (7 = 0) case, the A(k F ,u) has no 5-function at uj — 0. The 
spectral function goes as A(kp,uj) ~ (\u\ |ln(|o;|)| 3/2 ) _1 . In the dirty case 7 ^ 0, we have 
A(kp,uj) ~ I \n(\u\)\ 3 + 2ttZf5(uj). We can see by comparison of the 7 = and 7 = Ep lines 
in Figs. 2(c) that some of the spectral weight around the origin in the 7 = case has been 
transferred to the <5-function in the 7 7^ case. 

14 



In Fig. 01, we show our calculated Fermi distribution function [IT? 



n k = / —A{k,u) (2.22) 

J-oo 2,71 

for various values of the impurity scattering rate 7. We emphasize that 7 was included 
only in the dynamical screening function and not in the single-electron Green's function 
because we wanted to determine if the suppression of the emission of low energy virtual 
plasmons produces a discontinuity in n\.. In the next subsection, we discuss the effects of 
introducing disorder into the bare Green's function. Fig. f|(a) clearly shows a discontinuity 
in rik at k = kp for j/Ep 7^ 0. In Fig. f|(b), we show the calculated Zp as a function of the 
impurity scattering rate. For 7 = 0, Zp = indicating that there is no Fermi surface, but as 
scattering is increased Zp also increases until it saturates at very large 7 (where our results 
should not be trusted because our treatment ignores localization). Note that Zp goes to 
zero slowly as 7 — > 0, implying that even a small amount of impurity scattering results in a 
fairly pronounced discontinuity in n k at k F . In Fig. (|(c), we show the density of states (per 
spin) 

1 poo flh 

p(w) = sL^ ( * ,w) (2 - 23) 

for two different impurity scattering rates. [EJj] For a clean wire, there is a slow (inverse 



logarithmic) disappearance in p(oj) as uj — > 0, indicating the existence of a gap in the 
density of states at the Fermi surface. For non-zero 7, our numerical results indicate that 
p{uj) initially decays as uj — > 0, but then flattens out to a finite value, indicating the absence 
of a gap (in accordance with Fermi liquid theory). Furthermore, note that in this calculation, 
we have not included the effects of scattering in the Green's function in the self-energy (see 
subsection [11 D| ), which would smear the density of states further and also remove the gap 



associated with the singular behavior of a Luttinger liquid at the Fermi surface. 

Fig. ||(a) shows the inelastic scattering rates of quasiparticles in the conduction band 
T(k) = 2 |Im[E(&, u = £&)]! f° r parameters corresponding to a = 100 A and a density of 
n = 0.56 x 10 6 cm~ 1 in GaAs. Fig. |5](b) shows the corresponding inelastic mean free path, 
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I = v(k)/T(k), where v is the electron velocity. For 7 = 0, below a threshold wavevec- 
tor k c , there is no electron-electron scattering (within the RPA) because in a strictly one- 
dimensional system, conservation of energy and momentum restricts electron-electron scat- 
tering to an exchange of particles, which is not a randomizing process because electrons are 
indistinguishable. |?T| (Our treatment ignores multiparticle excitations, which will give rise 
to a nonzero scattering rate for k < k c .) For k > k c , a new scattering channel opens in 
which electrons genuinely emit plasmons (as opposed to the virtual plasmon excitations at 
the Fermi surface). The inelastic scattering rate diverges as ~ {k — k c )~ 1 / 2 as one approaches 
k c from above, due to the divergence in the density of states available for scattering right 
at the plasmon emission threshold, analogous to the divergence in the optical-scattering 



rate in one- dimensional systems. |g2| In section |V|, we discuss this singular behavior in the 
scattering rate in greater detail and propose that this, in principle, can be used to fabricate 
a novel device. For 7 7^ 0, the inelastic scattering rate remains finite because the plasmon 
line is broadened. Furthermore, the breaking of translational invariance relaxes momentum 
conservation, permitting inelastic scattering via single particle excitations for k < k c . 

In Fig. we show the results of the calculation of the the bandgap renormalization, 
which is the sum of Re[S(/c = 0, u> = £fc=o)] °f the conduction band electrons and the valence 
band holes. These results should be useful in explaining photoluminescence experiments 
in quantum wires, even though our calculation takes into account the screening effects of 
electrons only, while a photoexcited semiconductor in fact contains a finite density of both 
electrons and holes. This is justified because we expect the screening effect of the holes 
to have a negligible effect on the bandgap renormalization, since their large mass prevents 
them from screening effectively. 

Table [TI| summarizes the properties of a quasi one-dimensional quantum wire, for both 
clean and dirty systems. 
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D. Impurity effects on the bare Green function 



So far, we have ignored the effect of the impurity potential on the non-interacting Green 
function G (k,u); i.e., we have used the G (k,uj) for the impurity-free case to calculate the 
many-body properties of impure, disordered quantum wires. In fact, the disorder caused by 
the impurity potential in one-dimensional systems has a very significant effect, tending to 
localize all the eigenstates, P3,||| thus significantly changing Go(k, to) from the impurity-free 
case. In this subsection, we address the question whether using the non-disordered Go(k,u>) 
for the disordered case is justified. We also calculate G (k, uj) for the disordered case for 
some simple models. 

For a particular impurity configuration, assuming the exact non-interacting electron 
eigenstates ipi(x) with eigenenergies £j (with respect to the chemical potential), the non- 
interacting Green's function would be diagonal in the basis of these eigenstates, G (ij, to) = 
5ij{uo — £j + iO + sign(£j)) . Thus, for the non-interacting case, the distribution function at 
T = would be arii = 0(— i. e., the momentum index k in the pure case is simply replaced 
by the indices enumerating the exact eigenstates. Hence, our calculation corresponds to the 
situation where the k labels in the self-energy correspond not exactly to the momentum, 
but to the indices of the exact eigenstates of the disordered noninteracting system. Thus, 
the discontinuity at the Fermi surface is not a discontinuity in the distribution function 
in momentum space, but in the representation of the exact eigenstates of the disordered 
system, and our results should be interpreted in the basis of the disordered eigenstates. The 
arguments regarding the suppression of plasmons in the disordered system should still hold, 
since the fact that there are no low energy collective oscillations in disordered systems is 
independent of whatever basis is used. 

If one wanted to calculate Y,(k,uj) in the basis of the momentum eigenstates then the 
bare Green's function would have to include the effects of the impurity potential. In this 
subsection, we calculate the distribution function n^, which is the occupation of the disorder- 
free momentum eigenstates, for varying amounts of disorder, using two different methods. In 
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the first, we use the standard impurity averaged scheme p0|j45[| to calculate, to lowest order 



in the impurity scattering potential and in the absence of interactions between the electrons, 
the Green's function and the momentum distribution function n^. The Feynman diagram for 
this approximation is shown in Fig. [7[ We show that, as in higher dimensions, the presence of 
impurities broadens because the exact eigenstates are superpositions of different fc-states. 
Since there are no interactions between the electrons, the electrons still form a Fermi liquid, 
in the sense that there is a sharp discontinuity between the occupation of states above and 
below the Fermi energy. The second method we use to calculate n k involves putting in a 
term representing impurity scattering in the expression for the exact Green function of the 
Luttinger liquid. We show that the impurity scattering tends to remove the singularities 
associated with the non- Fermi liquid nature of the Luttinger liquid. A comparison of both 
momentum distribution functions indicates that for large enough impurity scattering rate, 
the noninteracting and interacting distribution functions are essentially indistinguishable, 
which implies that any singularities in n k associated with the interaction is masked by 
impurity effects. 

For the standard impurity self-energy for impurity scattering potentials of the form 
U 6(r - r { ), 

NAUnl 2 r°° dk 



S(fc, Uj) 



27T J-oo UJ — £fc — i0+ 

{|u; + uh 1/2 , ifcu + /i<0: 
(2.24) 
i{uj + Ai)" 1/2 , if u + fi > 0; 

where C = NJJq (j^j 2 = 7i m pV r ^F/2, where 7i mp is the impurity Born approximation scat- 
tering rate for an electron at the Fermi surface. Defining uj' = uj + /z, the spectral function, 
from Eq. ( ^Wj ) and ( ggf , is 



All i \ 27T C/y/u? 

A(k, u, = «, - Ai) = 1+ c M 3/ 2 ^ - ^) + „ V _ Ek? + C - 



(2.25) 



2 

where io§ is the solution of 



co - E k = j=. (2.26) 



The first term on the right hand side of the above equation is from a 5-function in the 
spectral function below the band-edge, corresponding to a bound state. 
The momentum distribution function, from Eq. ( |2.22| ), is 

- 1 du/_ CV^ 

Uk ~ 1 + CVo| 3 / 2 /2 + Jo vr uj'{uj' - E k f + C r { ' } 

The chemical potential /i is calculated self-consistently so that, as 7i mp is changed, the total 
density n = J dk n^j-a is kept constant. Fig. || shows /i as a function of 7i mp , and the 
corresponding momentum distribution functions. Note that, by keeping only the first order 
impurity scattering diagram, we are assuming weak scattering. Therefore, the results are 
unreliable in the large 7 imp (i.e., large disorder) limit, multiple scattering events are ignored. 
For example, if we included all the diagrams corresponding to repeated scattering off a single 
impurity, each such diagram with n impurity lines for a short-ranged interaction contributes 



a term NiUo(iUo^m/ (2uo + /i)) n , and therefore the entire series can be summed up to give 

NJI 2 

S(fc, u) = . (2.28) 

—iJ2uo/m — Uq 

Therefore for large Uq, the S(/c, u) is modified significantly from its lowest order form. There 
will also be significant contributions from scattering events off two or more impurities. 

We also investigate how disorder affects the distribution function n k of the Tomonaga- 
Luttinger model, where the kinetic energy of the electrons is given by the completely linear 
dispersion = vpk. The exact real-space Green function for the interacting Luttinger 
model where the interaction term is 



1 ^-^ TTVf + i 

H e-e = ~ }^ 9 — a k a k' a k'-q a k+g (2-29) 
q,k,k' " 



and where the interaction g cuts off for momenta larger that A (which would physically 
correspond to the band- width cutoff in a real system, and which we set equal to hp) can be 
found exactly, [|l7],ill and it is of the form 



G(x, t) = G (x, t)F(x, t), (2.30) 
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where F(x,t) is dependent on the interaction strength and bandwidth cutoff. 

The standard method for introducing impurities into the noninteracting Green's function 
is to introduce an electron lifetime r to the Green's function, 

G(k, u) = * . (2.31) 

to - 4fc + ^:sign(u;) 

which gives a real-space representation of G(x,t) of 



1 1 II 

GtapCM) = o ,.- n+ - ^ ^("9 )• (2 ' 32) 

27r s — i^i + zO + sign(t) v 2f r / 

Therefore, for the Green function in the case where both interaction and disorder exist, we 

make the plausible replacement of Go(x,t) by Gi mp (x,t) in Eq. (|2.30| ), yielding G(x,t) = 

Gi mp (x,t)F(x,t). From this, the momentum distribution function is computed to be 

1 1 f°° , sm(kx) \x\ . / , o, , . . 
n k = --— dx^-^expi-^-) (l + AV - a , 2.33 

2 27r j-oo x 2t>^r v 

where 

a = ((1+S) 1/2 - 1) 2 /(8(1 + ^) 1/2 ). (2.34) 

For g < 1, a w g 2 /32. 

The exact theories assume interactions which are constant (in momentum space) and 
short-ranged (in coordinate-space), which is clearly not the case for the bare Coulomb inter- 
action in a quasi-one-dimensional quantum wire. However, in many experimental situations, 
the quantum wires are surrounded by mobile charges from objects such as metallic gates 
or adjacent quantum wires which would tend to screen the Coulomb interaction between 
electrons within the same wire. This screening would reduce the range of the interactions, 
resulting in an effective intrawell interaction V e g(q) which is finite as q — > 0. The effective 
interaction for interaction for a wire screened by mobile carriers in an adjacent wire, for 
experimental parameters given in references 0, is on the order of g = 7.5, which gives 
a = 0.15 (see Appendix E). 

In Fig. [5] we show the the distribution function given by Eq. Q2.33p , for a = 0.15, 



and various values of 7. For 7 = 0, the distribution function has an inflection point, 
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I n& — 1/2| ~ \k — kp\ 2a a,s k kp, which takes the place of the discontinuity of the non- 
interacting Fermi liquid. The introduction of disorder smears this inflection point out, as it 
would with a Fermi surface in higher dimensions, and with increasing disorder it becomes 
indistinguishable from the shown in Fig. [8] for a non-interacting disordered system. Thus, 
impurity effects tend to mask the singular behavior of a Luttinger liquid, which would tend 
to complicate any experimental attempt to measure the the Luttinger liquid singularity in 
the nfc. 

E. Discussion of the zero-temperature self-energy 

In this subsection, we compare the relative merits of the Tomonaga-Luttinger model 
and our approach, and we discuss the justification for treating electrons in one-dimensional 
semiconductor quantum wires as Fermi liquids. 

As mentioned previously, interacting one- dimensional systems have been shown to be 
non- Fermi liquids, through the solution of the Tomonaga-Luttinger model. With the model's 
assumptions of infinite density of negative energy electrons in a completely linear dispersion 
and short-ranged interactions, all the vertex corrections in the self-energy are included and 
the solution is exact. On the other hand we assume a finite density of electrons in a parabolic 
energy dispersion and we use the actual Coulomb interaction |2(J between electrons for a 
rectangular well, but we carry out only the leading order self-energy calculation in the dy- 
namically screened interaction. Thus, previous work has concentrated on the exact solution 
of an artificial model, whereas our work is an approximate solution of a realistic model. 

While the solution of the Tomonaga-Luttinger model is mathematically exact, it is dif- 
ficult to see from the exact solution which physical processes are responsible for the drastic 
changes in the low-energy behavior of the system and the disappearance of the Fermi sur- 
face. Furthermore, since the model linearizes the spectrum about the Fermi surface, it is 
not reliable in describing physics at energies far away from the Fermi energy, when the band 
parabolicity becomes important. While our calculation is only approximate, we can use it 
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to isolate the processes that have a a decisive role in determining the low-energy properties 
of the one-dimensional interacting electron gas. We have shown that, physically, the virtual 
plasmon emission at low energies is responsible for the disappearance of the Fermi surface. 
In Appendix F, we show that inclusion of leading order vertex corrections in the form used 
by Mahan and Sernelius [071 does not change our results at the Fermi surface. We can also 



explore the high energy processes, such as inelastic scattering of hot electrons, because we 
have included band parabolicity and we can calculate the band-gap renormalization effects, 
because unlike the Tomonaga-Luttinger model, our model has a band minimum. 

We conclude this section by discussing the viability of treating one-dimensional elec- 
tron systems as Fermi liquids. One-dimensional electron systems in solids have theoretically 
strikingly different behavior from their higher dimensional counterparts. In addition to the 



Luttinger liquid behavior, there are other effects, such as the Anderson |f43| , f44|| localization 
in the presence of disorder and the Peierls lattice distortion in the presence of electron- 
phonon coupling, which theoretically should drive a one-dimensional electron system away 
from Fermi liquid behavior. It would seem to be completely inappropriate to think about 
one-dimensional electron systems as Fermi liquids. We argue, however, that in real semi- 
conductor quantum wires, various factors may serve to stabilize Fermi-liquid behavior. 

In one dimension, for non-interacting electrons in the presence of any disorder (again, 
invariably present in real wires), all states are localized and the concept of a one-dimensional 
electron gas does not (except in the unrealistic no-impurity idealization) apply. Therefore, in 
the absence of electron-electron interactions, all the quantum wire electronic states are ex- 
ponentially Anderson-localized. The situation is again fundamentally different from higher 
dimensional electron systems where finite disorder strength is required to strongly localize 
the electron states. We argue, however, that in the state of the art semiconductor quantum 
wires the effect of Anderson localization is negligibly small because typically the localization 
lengths are very long (many microns) due to the high quality of the sample and the mod- 
ulation doping technique. Furthermore, the effects of electron-electron interaction on the 
localized states of a disordered one- dimensional system are by no means completely under- 
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stood, and theories exist which indicate that small amounts of disorder may not localize an 



interacting one-dimensional electron gas. |E8] In any case, for weakly disordered quantum 



wires, the screening effect of the electrons would reduce the effect of the impurity potential 
and increase the localization length. Thus, so long as the localization length is much longer 
than the length scale of the experimental probe (e.g., in Raman scattering, that would be 
the wavelength of the photons), the electrons in the quantum wires may be considered to 
be extended for all practical purposes. 

The presence of coupling of electrons to the acoustic-phonons produces, in the ideal zero- 



temperature situation, a Peierls distortion, [49 in which the lattice distorts to lower the 



energy of the electron gas at the expense of a smaller increase in the lattice potential energy. 
The Peierls distorted state has a bandgap at the Fermi wavevector, and so theoretically quasi- 
one-dimensional semiconductor quantum wires should be insulators at T = 0. However, 
the electron-phonon interaction via the deformation potential coupling is so small that the 
Peierls transition temperature is on the order of 10~ 3 K, or much less, depending on the 
density. []5D| Furthermore, the transition is suppressed altogether if the impurity scattering 



rate is larger than the transition temperature, and therefore in semiconductor quantum 
wires, since the Peierls transition temperature is so low, even a small amount of disorder 
will remove the transition. 

Finally, we have shown that the impurities serve to remove the mechanism (virtual 
plasmon emission) responsible for destroying the Fermi surface. Therefore, as with the 
Peierls distortion, the presence of impurities seems to prevent the system from evolving away 
from Fermi liquid behavior. Thus, for all practical purposes, real semiconductor quantum 
wires can be considered to be normal one- dimensional Fermi liquid systems. 

III. MANY-BODY PROPERTIES OF QUANTUM WIRES 
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AT FINITE TEMPERATURE 



In general, many-body calculations in metals [52] have used the zero-temperature for- 
malism because the energy scales intrinsic to the problem (Fermi energy, plasmon energy, 
etc) are usually much larger than the temperature T. In contrast, in semiconductors, espe- 
cially in artificial structures of reduced dimensionality, because of the low electron densities 
and large lattice dielectric constants involved, the experimental temperature can be com- 
parable to the intrinsic energy scales of the electron gas. For example, for a quantum wire 
density of 5 x 10 5 cm _1 , the Fermi energies are of order 5 meV ~ 50 K. Furthermore, the 
one-dimensional plasmons are gapless and have a dispersion that is almost linear, implying 
that there is a large density of low-energy collective excitations. Thus, even at relatively 
low temperatures, the intrinsic energy scales of the electrons and of collective excitations 
can be comparable to the temperature. Therefore, the zero-temperature formalism may 
not provide an adequate description of the system, and the finite-temperature formalism is 
needed. |kJ In this section, we calculate the finite temperature self-energy, again using the 
GW approximation (i.e., the Feynman diagram shown in Fig. 0(a)). 

In the T = case, by obeying the Feynman rules, one obtains an integral expression for 
H(k,u) in the GW approximation. To numerically evaluate to), we used the so-called 
"line and pole" decomposition, which is obtained by a contour deformation of the original 
integral expression for u). [[18 , 27] The "line and pole" decomposition is more efficient for 
the purpose of numerical computation than the original expression. However, to the best of 
our knowledge, this "line and pole" decomposition has not been generalized, in the published 
literature, to finite-temperature self-energies. In this section, we introduce a general method 
for obtaining expressions for the real-frequency finite-temperature self-energies, based on 
analytic continuation of the temperature (Matsubara) Green's function, which, in the case 
of the GW approximation, yields the finite-temperature generalization of the "line and pole" 
decomposition. This method is valid for arbitrary dimensionality. We present the results of 
our calculations for a semiconductor quantum wire using this formalism. 
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A. Formalism 



Matsubara [|53J first introduced a many-body formalism for finite-temperature systems, 
which is formally identical to the many-body zero-temperature formalism. While the finite- 
temperature formalism is easier to handle than the zero-temperature formalism in some 
ways, there is one added complication. With the finite-temperature formalism, one obtains 
an expression for the electron self-energy a(iu n ) that is only valid at discrete points on the 
complex frequency plane, iv n = i(2n + 1)ttT, where n is an integer. However, the quantity 
that is of interest for physical properties of a system is the self-energy at real frequencies. 
Therefore, the u{iv n ) must be analytically continued to the complex plane to obtain the 
self-energy that is valid for all complex frequencies, |34| (for convenience, at this stage, 
we have suppressed all arguments of the self-energy except for frequency) from which the 
real-frequency (retarded) self-energy, the quantity relevant to experiments, can be obtained 
by setting z = u + i0 + . Note that in this section, we use the convention that functions 
denoted by upper case characters are analytic in the frequency variable, while those denoted 
by lower case characters may be nonanalytic. 

In principle, from a(iv n ), one can obtain a formal expression for £(z) in terms of integrals 



over spectral representations. [18] For example, in the GW approximation, the expression 



for the real-frequency retarded self-energy is 

E ( k, w) = /Ar ^ - B Y } . n+ (M^)+MM (3-D 

where al is the dimensionality of the system, B(q,u) is the spectral function of the screened 
Coulomb interaction, given by 

B{q, oj) = 2Re[W(q, u)] = -i(W(q, u + i0 + ) - W{q, to - i0 + )) , (3.2) 

and 

UF{X) = exp^/T) + 1 (3 - 3a) 

nB{x) = exp(x/ 1 T) - 1 (3 ' 3b) 
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are the bose and fermi functions. Expressions for £(o>) obtained in this manner, however, 
involve integration over one or more frequency variables, which makes them inefficient for 
use in numerical computations. On previous occasions various approximations such as the 



plasmon-pole approximation were used to obtain the finite-temperature self-energy. |j3] The 
plasmon-pole approximation assumes that the spectral function of the screened Coulomb 
interaction, B(q,u), can be replaced by a ^-function, whose weight and position are given 
by sum rules. The ^-function removes the frequency integration, considerably simplifying 
the numerical computation. Of course, this is an uncontrolled approximation, and it has 
been shown that on some occasions it is quantitatively inaccurate, [|55[] and it is desirable 
to evaluate the expression Eq. ( j3.1|) efficiently, without making any approximations. In this 
paper, we present a more direct method of analytically continuing the a(iv n ) to £(z) which 
yields an exact expression which is much more amenable to numerical calculation than the 
integrals over spectral functions. We then apply this method to calculate various many-body 
properties of electrons in a quantum wire. 

This method is based on the principle that the the analytic continuation of a(iv n ) satisfies 
certain necessary and sufficient conditions, and once these conditions are satisfied, one is 
ensured that one has obtained the unique analytic continuation, E(,z). These conditions 
are as follows. (1) £(z) is analytic on the entire complex frequency plane, with the exception 
of possible branch cuts on the real axis @] (henceforth, when we say a function is "analytic" 
it is with the implicit understanding that it could have branch cuts on the real axis); (2) 
E(z = iv n ) = a{iu n ) for all iv n \ and (3) £(z) goes to a constant as \z\ — > oo. The method is 
based on systematically fulfilling each of the above conditions one by one, leading us directly 
from a(iu n ) to its analytic continuation £(z). 

To expound this method, we discuss how it can be applied to a simple example, that of 
the self-energy within the GW approximation (see Fig. 0(a)) to obtain the finite-temperature 
generalization of the "line and pole" decomposition. From the Feynman rules, the expression 
for the Matsubara self-energy for imaginary frequencies is [ITS] 
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a{k,iu n ) = - / j—rj T^——, — -- — , (3.4) 

where the summation is over the boson frequencies iu n = i2irnT (n are integers). This 
expression can be written as a sum of a frequency-independent (and hence analytic with 
respect to frequency) exchange and a frequency-dependent correlation part, fl"8| 



where 



cr(k, iu n ) = cr cx (k) + cr CO r(k, iu n ); (3.5a) 

^ex(k) = -J ^ K(q) (6*k)> (3.5b) 

/dq 

h^iun), (3.5c) 



tZ lV n + M n - £ k +q 



Here, 



w(q, lUJ n ) 



^e(q, iu n ) 

is the difference between the screened and bare Coulomb interactions. The problem is 
to analytically continue <7 cor (k, ii/ n ) to S cor (k, z). The most obvious attempt, the simple 
substitution iv n — > z in Eq. Q3.6|) , does not give the desired analytic continuation because 
there are poles in the integrand h^{z) at z — £k+q — iwn> which translate into branch cuts 
in cx cor (k, z) at lm[z] = 27rnT for all n. This clearly violates the condition (1) stated above. 

The cr cor (k, iv n ) is obtained by an integral over q of an imaginary-frequency-dependent 
integrand h^ ^ii/n), as shown in Eq. ( |3.5c|) . Thus, if we obtain analytic continuation H^ ^z) 



of h^qfon), then 



E OT (M = -/ j^H^z), (3.8) 



automatically satisfies conditions (1) and (2) above for the analytic continuation of H(k,z). 
f57| Hence, our problem is reduced to constructing function H^q(z) such that (1) it is 
analytic (in the sense mentioned above, i.e., analytic off the real frequency axis), and (2) 
Hk,q(z = iv n ) = h k> q(iv n ). 
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The outline of the procedure for obtaining the function H^q(z) which satisfies conditions 
(1) and (2) is as follows. First, we write down a function H£ (z) = hk i(l (z) + hk M (z), where 
hk,q(z) is chosen so that it cancels all the singularities in hk iCl (z) on the complex plane. 
Below, we describe a systematic method of obtaining H£ (z) from h^^z). Thus, H£ (z) 
is analytic, fulfilling condition (2). However, hy_ ;Cl (z) is in general nonzero at z — iv n , and 
hence H^^iv^) ^ ft<k,q(^n)j violating condition (2). The second step is therefore to add 
an additional analytic term H^^z) which cancels h^^z) at all z = iu n . Since the function 
Hk,q( z ) + H^ ^z) is analytic and equals hk )Cl (z) for all z = iu n , fulfilling both conditions 
(1) and (2), it is the desired analytic continuation H^{z). Condition (3) can be checked 
in the end; in the case of the GW approximation (and in other cases we have studied) it is 
satisfied. 

In the case of the GW self-energy, the expression obtained from the Feynman rules is 

hka(z)=T}^—— . (3.9 

Z + lU n - £k+q 

One can construct an analytic function which contains hk M (z) as one of its terms by writing 
a contour integral 



duj w(q, uj)n B {uj) 
c 2ni z + oj — £ k+q 



<(*)= /,S / r , (s-io) 



where the contour of integration C is shown in Fig. |H| H£ (z) is clearly analytic (off the 
real axis) in the variable z, since the function (z + u — £k+q)~ x has no poles off the real axis. 
By the residue theorem, the integral Eq. ( |3.9| ) can be evaluated in terms of the residues of 
the poles contained within the contour C. This gives H£ (z) = h^{z) + h^{z). The term 
^k,q(^)) comes from the poles of as can be seen more clearly by writing tib(w) as the 

sum of its poles, 

M0 = 4 + r£-r^rr- (3.ii) 

2 C — iu„ 

In addition to the poles of the bose function, there is also a pole due to the denominator of 
Eq. ( |3.9|) which gives the second term 

28 



^k,qO) = w(q, ^k+ q - z) n B (£ k+q - z). 



(3.12) 



(Note that w(q, u;) is analytic everywhere except for a branch cut on the real axis, and so it 
does not have poles which contribute to H^ q (z).) Despite being the sum of two nonanalytic 
functions, H^Jz) is analytic because the singularities present in both the functions precisely 



cancel each other. This can be seen clearly if we use Eq. ( p. lip to expand the bose function 
of hu,q(z), yielding 



•-w(q, Ck+q - z) + T ^ — ^ 

2 Z + lU n - £ k+q 



(3.13) 



which explicitly shows that both the denominator and numerator vanish at z = £k+q — ioj n . 

However, because h ktqi (iu n ) ^ 0, if k+q (iz/ n ) ^ ^k,q(^n) an d hence H^^z) fails to fulfil 
condition (2). This can be remedied by adding an analytic function which cancels /i k , q (z) 
at all z = iv n . Noting that %v n = i(2n + 1)ttT, for all integers n, 



n B (^+<i - iv n ) = exp(£ k+q /T) exp(-(2n + 1)tt) - 1 



exp(£k+ q /r) + 1 = -n F (£ k+q 



(3.14) 



and thus h ktti {iv n ) = — w(q, £k+ q — iv n ) n F{Ck+q), implying that the analytic term needed 
to cancel h k! q(iu n ) is 



H L M (z) = w(q,£ k+q - z)n F (£ k+q ). 



(3.15) 



Hence, H kM (z) = H^^z) + iJ k q (z) is given by Eqs. (|3.13|) and (|3.15| ), and the correlation 
self-energy in the GW approximation is 



^cor(k, Z) — 



z + iu n - & 



k+q 



/ (2^F ^ k+q ~ ^ ( ni? ^ k+q ~~ z ) + 



(27T) 



(3.16) 
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The retarded self-energy is obtained by setting z — > u+iO + . The first and second terms on 
the right hand side of the last equality of Eq. ( p,16| ) are, respectively, the finite-temperature 
generalization of the so-called "line" and "pole" components of the GW approximation of 
the T = correlation self-energy. [^] As in the zero-temperature case, the line contribution 
is completely real because u>(q, — iu> n ) = w*(q,iu n ) and hence the total contribution to the 
imaginary part of E(k, u) comes from the pole part. 

As in the T = case, in the GW approximation the "on- shell" imaginary part of the 
self-energy, |Im[£(k, to = £k)]|, is na ff the sum of the Born-approximation electron and hole 
scattering rates. The Born-approximation rate for a particle to be scattered with change in 
momentum q and energy E in the particle is |58 



P(q,E) = -2V c (q)lm[e-\q,E)}n B (E) 



(3.17) 



Using the identity 



ns(6c+q - k>) + ^F(6c+q) = ™I?(£k+q ~ Oj) \ - ^(Ck+q) ~ U B (uJ - £k+q)™F (£k+q) ) (3-18) 



we can write 2 |Im[£(k, £k)]| = 7 e (k) + 7h(k) where 

cfq 



7e(k) 



7h(k) 



(2ix) d Im ^ ^ k+q _ ^ nB ^ k+q _ ^ i 1 ~ nF (&+i 

c?q r 

d-P(q,fk+ q ) 1 -n F (£ k+q 



(27T) 

(2n) d 
(2n) d 



y c (q)Im[e x (q, £ k+q - 6c)]n B (6c - fk+qWCfk+q) 

^(q,^k - ^k+q)^F(^k+q) 



(3.19a) 



(3.19b) 



are the total Born approximation electron and hole scattering rates, respectively. |5^] The 
factors of 1 — ni?(£k+q) and np(£k+ q ) that multiply the Born approximation scattering rates 
in Eq. ( |3.19b| ) are due to the Pauli exclusion principle. An electron (hole) can scatter to a 
final state k + q only if it is unoccupied (occupied), and the factor 1 — np(£ k + q ) (np(£k+ q )) is 
the probability that the final state is unoccupied (occupied) by an electron from the system. 

The method outlined above which we have applied to the GW approximation can also be 
applied to higher-order diagrams, although the procedure quickly becomes considerably more 



30 



complicated. In a higher-order self-energy with an external Matsubara Fermi frequencies iu m , 
the Feynman rules result in an expression with a summation of a function f(iu m , u ni , u nN ), 
over several boson frequencies u ni , ...,u nN , 

h{iv m ) = ■■■ f{iv m ,u ni ,...,uj nN ) (3.20) 

In f(iu m ,cu ni , ...,u nN ), there are several bare Green's function which contain both %v m and 
one (or more) internal Bose frequency iuj ni ] e.g. terms of the form (iuj n .-\-iv m — ^k+q) -1 - (One 
can always write down the self-energy in a form such that the external frequency occurs only 
in the bare Green's functions, and not in the screened Coulomb interactions.) To obtain the 
analytic function H A (z) which contains h(z) as one of its term, one writes down an integral 



H (z) = dux... du N n B (u 1 )...n B (uj N )f(z,uj 1 ,...,uj N ), (3.21) 

J Ci JCn 

where every summation over iou ni , has been replaced replaced by an integral over u>i, iv m is 
replaced by z and / is multiplied by the a product of bose factors. Since z only appears in the 
integrand in the form of factors (u + z — £k+q) _1 > H A (z) is analytic off the real axis. By using 
the residue theorem, one can evaluate the integral in terms the poles of the integrand. The 
nsi^i) terms give summations over Bose functions, which correspond exactly to the original 
expression, Eq. ( |3.20| ), with iv m replaced by z. There are also other terms generated which 
serve to cancel the singularities of Eq. ( p.20| ) on the complex plane. As in the example cited 
above for the GW approximation, this function fails condition (2), and therefore analytic 
functions must be found that cancel the unwanted terms on the points z = iv m . Generally, 
this can be done by inspection of the resulting term, and by use of Eq. ( |3.14| ). We have 
used this method to obtain an expression for the self-energy which is second-order in the 
screened interaction, which we give in Appendix G. 

Alternatively, one can also derive the finite-temperature "line and pole" expression of 
self-energy from the spectral representation of the GW approximation of the self-energy, 
given by Eq. Q3.1| ). The integral over the spectral function for the Coulomb interaction can 



be rewritten as a contour integration over the contour indicated in Fig. [10], with the spectral 
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function in the integrand replaced by i~ l w(q, to). Then, by the residue theorem, one obtains 
precisely the expression obtained above in Eq. ( |3.16| ). This procedure also works for higher 
order diagrams. One first obtains the expression for the self-energy in terms of integrals over 
spectral functions, and then by writing these integration of spectral functions as contours 
of the form shown in Fig. |TU] and evaluating the contributions due to the poles within the 
contour, one obtains an expression involving sums over imaginary frequencies. 

The main advantage of the expressions generated by the method outlined above is that 
the frequency integrals in the spectral representation method are generally replaced by 
imaginary frequency summations, which are numerically more efficient to perform, especially 
at high temperatures. Unfortunately, with the higher order diagrams, there are terms that 
still involve frequency integrations, but the number of these frequency integrations is reduced 
from the spectral representation expression of the self-energy. 

B. Results for quantum wire 

We applied the formalism developed above to electrons in a semiconductor quantum wire. 
As mentioned previously, the densities of current quantum wires are necessarily low, and the 
Fermi energies are on the order of 5 meV ~ 50 K. Thus, even small temperatures may affect 
their many-body properties significantly. We applied Eq. (|3.16|) to the calculation of the 
self-energy, spectral function and band-gap renormalization of electrons in a one-dimensional 
quantum wire in the extreme quantum limit [i.e., assuming that the electrons only occupy 
the lowest energy subband). 

We used the RPA form for the dielectric function e(q,z) = 1 — V c (q)U (q, z). Using the 
familiar expression for Uo(q,z), given by Lindhard, |60|JT8|] to numerically compute the real 
part of the polarizability would have been complicated by a principal part divergence in the 
integrand. Therefore, instead, we used 



32 



U.o(q,z;T,n) 



Je-^IT) (x + 1 V 



-wn (x + If 

7o (x + l) 2 (3_ 22 ) 



/• e<M/T) 

/ 7 — —r^ Tl o(q,z;T = 0,fi-T\n(x)), if // < 0; 



which is a computationally efficient version of the expression given by Maldague |6TJ] for 
the finite-temperature RPA polarizability. In performing the summations over imaginary 
frequencies of the "line" part of Y>(k,uj), we used the fact that the asymptotic behavior of 
the components go as n~ 5 / 2 as n — ► oo (the asymptotic form for these components for an 
infinitesimally thin d = 2 electron gas is n~ 2 and for d — 3 is n~ 5 / 2 ), which allowed us to 
use a slightly modified form of the Richardson extrapolation procedure |62| to obtain rapid 
convergence of the sum. Also, since there are mutually canceling principal part divergences 
(in the g-integration) in the n = term of the "line" part and the real part of the "pole" 
term, it is advantageous to numerically evaluate both these terms together. 

We calculated the real and imaginary parts of the self-energy and the spectral function 
of a quantum wire as a function of frequency, for several temperatures. The density was 
kept constant by adjusting the non- interacting electron gas chemical potential /io(^) so that 



the integral / dk n-p(Ek — Ho(T))/tt was kept constant. We show our results in Fig. |TT 
The discontinuities and kinks in Im[E(a>)] at T = 0, which arise from virtual plasmon and 
single-particle emission thresholds, broaden with increasing temperature because the plas- 
mon peaks are broadened by Landau damping, and the boundaries of the single-particle 
continuum are thermally smeared. At non-zero temperatures, a logarithmic divergence de- 
velops in the imaginary part of the self-energy at u = Ek — fJ,o(T)], and is accompanied 
by a discontinuity in Re[E(fc, £&)], since the real and imaginary parts of S(fc, u) are related by 
the Kramers-Kronig relations. This divergence is due to the non-integrable l^l -1 singularity 
as q — > in the integrand Eq. ( |3.16| ), which arises from the product of the bose factor, which 
goes as g -1 , and the Im[e _1 (g,^ +g — which goes as sign(g). Physically, this singularity 
corresponds to the singular divergence in the Born-approximation electron-electron scatter- 
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ing at small momentum transfer in d = 1, and is unique to one- dimensional systems, since 
in higher dimensions, the singularity in the integral is removed by the phase-space factor of 
q d ~ x . It is not obvious whether this singularity exists in the exact u) or is an artifact 
of the approximations we have used. 

In order to calculate the spectral function, one needs to know the finite temperature 
chemical potential for the interacting system. At T = 0, the chemical potential /i of the 
interacting system is simply given by Re[E(fc j p, uo = 0)] (as is required to form a discontinuity 
or a singularity at the Fermi surface). Unfortunately, at finite temperatures, there is no such 
obvious prescription to find fi(T) for the interacting system. In principle, one could search for 
the /x(T) which allows the sum rule Eq. ( |2 . 1 7| ) to be satisfied for all k (we expect Eq. ( [2.171) 
to be satisfied because we are using consistent particle-conserving approximations |kJ). 
However, in practice, the integrated spectral weights are not particularly sensitive to changes 
in the chemical potential and we were not able to determine \x accurately in the attempts 
we made. By the same token, the spectral functions themselves are not very sensitive to 
the choice of chemical potential. We used the T = value of the chemical potential of the 
interacting system, with a temperature correction based on the temperature dependence 
of iiq(T) for the non-interacting system. As expected, the spectral functions broaden with 
increasing temperature, corresponding to an increase in the quasiparticle decay rates caused 
by the presence of thermal excitations. The dips at uj = which go as ~ | ln(u; — ^fc)| _1 , 
which are brought about by the logarithmic divergences of Im[£(A;,a>)]. 



In Fig. [12], we show the electron and hole band-gap renormalization due to the presence of 
conduction electrons alone. We note again that holes are not expected to change the results 
significantly because holes screen weakly. Due to the discontinuity in Re[£(fc, we take 
| Re[£(fc,£/ c + + )-|-X!(fc,£fc — + )] at k = to be the band-gap renormalization. We find that 
for very low densities, where the Fermi temperature is low, the band-gap renormalization can 
change by approximately an order of magnitude when the temperature increases from T = 
K to T = 300 K. Such changes should be measurable in photoluminescence experiments. 
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IV. DEVICE POSSIBILITIES FOR QUANTUM WIRES 



Part of the reason why so much excitement has been generated by quantum wires is the 
possibility that they might have novel device applications. In the early eighties, Sakaki |T4 



argued that electrons in quantum wires could have extremely large mobilities because elastic 
collisions with small momentum transfer are suppressed in one-dimensional systems, which 
would make them very attractive for use in semiconductor devices. More recently, electronic 
devices based on the quantum nature of electrons and on the single-particle Coulomb 
charging energies in mesoscopic systems p5J have been proposed. [BH] In addition, various 



proposals for uses of quantum wires as lasers and optical switches have been advanced. |67 



In this section, we propose another interesting device principle of a one- dimensional quantum 
wire, based on the the many-body properties which are peculiar to quasi-one-dimensional 
systems. 

The device principle we propose is based on the many-body properties an interacting one- 
dimensional Fermi system. We show that it may be possible to obtain a device with large 
and sudden onset of negative differential resistance (NDR) (i.e., dl/dV < 0). This sudden 
onset of NDR could be exploited to produce a transistor, while the NDR itself suggests 
that this device might be used as an oscillator (e.g., in analogy with the Gunn oscillator 
or, more recently, the resonant tunneling diode In the proposed device principle, 



the predicted NDR is associated with a sharp change in the inelastic mean free path, due 
to electron-electron interactions (or more specifically, electron-plasmon scattering) of the 
injected electrons at a specific injection energy — in the ideal system at T = 0, the mean 
free path changes from being infinite below the threshold energy to being zero above it. 

The device principle which we propose may be experimentally observed in the quasi- 
one-dimensional version of the tunneling hot electron transistor amplifier (THETA), shown 
schematically in Fig. |H|(a), which has been fabricated successfully in three and two dimen- 
sions. JFOJ As in previous sections, we assume that the quasi-one-dimensional device is in 
the extreme quantum limit; i.e., that all the electrons are in the lowest energy subband 
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in the device. Electrons are injected from an emitter at energies above the Fermi energy 
Ep into a base region which contains (either through doping or electrostatic confinement) a 
one-dimensional electron gas, and the injected electrons that travel through the base region 
enter the collector on the opposite side of the base. The fraction of electrons a that reach 
the collector goes as a ~ e~ d /\ where d is the width of the base region and I is the electron 
mean free path. The I is equal to VkT(k), where Vk is the electron velocity and T(k) is the 
inelastic scattering rate, which is a strong function of the injection energy of the electrons 
and the electron density in the base region. Thus, by varying the scattering rate by changing 
the injection energy, one can achieve a significant change in the electron mean free path and 
hence the emitter-collector current. 

In two and three dimensions, the main scattering mechanism for these electrons in the 
THETA devices are the coupled plasmon-optic-phonon modes. [FT] However, in the semi- 



conductor quantum wires in the extreme quantum-limit that are currently being fabricated, 
the densities of the electrons in the base are so low that all the energy scales associated with 
the electron gas and operation of the device (Ep, plasmon energy and electron injection en- 
ergy) are much smaller than the optic-phonon energy, and therefore the optic-phonons play 
a negligible role. Acoustic phonons couple very weakly to electrons in III-V semiconductors, 
|72fl and the associated scattering rates are on the order of 10 10 s _1 and can be ignored when 



compared to the scattering mechanism discussed in this paper. We assume, for the purpose 
of this paper, that impurity scattering in the wires is negligible, which is not unreasonable 
given the excellent and continually improving techniques for fabricating these mesoscopic 
systems. This last assumption is equivalent to assuming that the elastic mean free paths 
are much longer than the inelastic mean free path to be calculated in this paper — given 
that our calculated inelastic mean free paths are generally a few thousand A or less, and 
in good quality quantum wires, elastic mean free paths are many microns, the neglect of 
impurity scattering is a good approximation for our purposes. Furthermore, the restriction 
of the scattering phase-space in one-dimension further reduces the scattering rates of both 
impurity and acoustic phonons. fl4| Thus the main scattering mechanism for an injected 
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electron is the interaction with the electron gas in the base. 

We use the Born approximation in our calculation of the scattering rate of the injected 
electron with the electron gas. The Born approximation calculation of the scattering rate is 
equivalent to the calculation of Im[E] in the GW approximation described in earlier sections 
of this paper. This approximation includes single-particle scattering and scattering with col- 
lective modes of the plasma, but ignores multi-particle excitations, which are expected to be 
smaller because they are higher order in the screened interaction. However, as we have men- 
tioned in the previous section, in one dimension and at finite temperature, a straightforward 
integration of the Born approximation scattering probability over all wavevectors yields a 
total Born scattering rate that is infinite, due to the presence of a non-integrable divergence 
in the number of small wavevector scattering events. These small-wavevector scattering 
events are associated with single-particle scatterings (in which the injected particle scatters 
with an individual electron in the electron gas in the base). The problem is that the Born 
approximation treats the injected particle as an external distinguishable particle, and gives 
the rate at which this particle is scattered out of its initial state. However, in actual fact, 
the injected electron is indistinguishable from the other electrons in the base, which is very 
significant in one-dimensional systems, because a single-particle electron-electron scattering 
event involves an exchange of electrons, implying that the final state after the scattering 
is exactly the same as the initial state. Therefore, a scattering rate associated with single- 
particle electron-electron collisions in one-dimension is spurious because these "collisions" 
do not change the state of the system. It is the transfer of momentum of the injected the 
particle to the plasmons of the system that degrades the current that flows from the emitter 
to the collector. The plasmons are collective modes associated with the electrons in the base 
and cannot exist outside the base region, and hence cannot carry current into the collector. 
Therefore, to determine the mean free path of the nonequilibrium distribution of electrons 
injected into the base, we should only include the scattering rate of the injected electron 
with plasmons. Unfortunately, within the Born approximation at finite temperatures, there 
is no way of separating the spurious single-particle scattering from the plasmon scattering, 
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since the sharp boundaries and line-widths of the single-particle continuum and plasmon 
dispersion which exist at T = both become thermally smeared and merge into each other 
as the temperature increases. ||75 ] 



Our solution to this problem is to calculate the momentum scattering rate, which is the 
integral of the Born approximation probabilities J57] weighted by the loss of momentum, 



T t (k) = 2 f° dq \ V c (q) Imle- 1 (q,uj k (q))] n B {u k {q)) [1 - n F (£ k+q )], (4.1) 

where k is the momentum of the injected electron, io k (q) = — £fc- g is the energy lost by the 
electron when is loses momentum q, and e(q,u) is the dielectric function calculated within 
the RPA. T t (/c) is generally the quantity that is relevant for transport, since it measures the 
degradation of the electron current. The form Eq. ( f4.1|) negates the effect of the spurious 
small q single-particle scattering events by giving them small weights and by canceling 
momentum loss against momentum gain. The main contribution to T t (k) should come from 
coherent peaks associated with the scattering of plasmons. In Fig. |TJ], we show the results 
of our calculation of T t (/c), and the corresponding inelastic mean free path, l k = v k /T t (k), 
for various temperatures. 

We concentrate on the zero-temperature T t (k), the simplest case. At T = 0, there 
are no thermally excited particles that contribute to the spurious single-particle scattering 
rate, and hence the only contribution to T t (/c) comes from the plasmon emission. However, 
not all injected electrons can emit plasmons. Because the plasmon dispersion in quasi-one- 
dimensional systems for small q goes as to(q) ~ q\ \og{qa)\ 1 ^ 2 , where a is the width of the 
wire, only injected electrons with large enough kinetic energies can emit plasmons (see Fig. 
|i~3|(b)). For a given density n, there is therefore a threshold wavevector k c (n) below which 
no plasmon emission can take place. Within the approximations we have used and at T = 0, 
as k is increased through k c , the scattering rate jumps from zero to infinity (equivalently, 
the mean free path falls from infinity to zero). The divergence in T t (k) as k — > k+, is 
(k — kc)~ l l 2 . Below, we derive this form of the divergence, and we show that under very 
special circumstances, this divergence may be stronger. 
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At T = 0, the electrons can only scatter by plasmon emission if momentum and energy 
conservation of the system is obeyed (and the final state of the electron k — q is initially 
unoccupied). One can graphically determine if it is possible to emit plasmons by plotting, 
on the same graph, the plasmon dispersion curve, uj p (q) [given by e(q,u p (q)) = 0], and 
the energy-loss vs. momentum- loss curve for an electron with initial momentum k, u^q) = 
£,k—£,k~q- Energy and momentum conservation of the system is obeyed and plasmon emission 
is allowed only when the curves u> p (q) and uJk{q) intersect. The wavevectors q* at which oJ p (q) 
and uik{q) intersect correspond to wavevectors of plasmons which the electrons are permitted 
to emit. The scattering rate at T = due to plasmon emission is given by 



del q, LUk(q) l - \ 

r P (fc) = EW(*)) { ^ — g l qHk y ( 4 - 2 ) 



dq 

where the summation is over all intersections q*. We are interested in finding the T = 

behavior of T p (k) near Jv til 6 minimum injected electron momentum at which the 

intersection between u k (q) and ui p {q) occurs. 

Defining q c and u c as the wavevector and energy of the at which the intersection of w p (q) 
and Uk{q) occurs for k = k c , the following relations hold: 



lo c = 

m 



-(kcq c - l -q 2 X (4.3a) 



e(q c ,uj c )=0, (4.3b) 
du kc (q) _ du p (q) 



(4.3c) 

<?=<?c 



dq i=<ic dq 

The first two equations are immediate consequences of the above definitions. The last 
equality, Eq. (|4.3cj ), results from the fact that, at k = k c , when uik{q) just impinges on the 



CG> p (g), the slopes of the two curves are equal (as can be seen by an inspection of Fig. |l"3"l(b)). 
Let co>i, qi and k± be the deviations away from uo c , q c and k c , 

q = q c + qu (4.4) 
k — ~\~ k \ . 
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In terms of q\ and k\, the energy vs. momentum loss curve is 

u k {q) =uj c + —(qoh + (h - q )q 1 + k x q x - -q(). (4.5) 

The dielectric function can be expanded around q c and u c in the form 

e(q, to) « -aq 1 + bu x + 0(q\, u\, q\U\), (4.6) 

where a, b > (this inequality can be deduced from the asymptotic behavior of e(q, &)). Eq. 
( [4.6| ) implies that the slope is at uj p (q c ) is a/6 and hence, by Eq. ( [4.3cj ), 

m b 

Substituting the expression for Uk(q) given by Eq. ( f4.5| ) into Eq. ( [4.6| ) yields 

e(q, UJk{q)) = -aqi + — Iqoh + (k - q )q 1 ] + 0(qf, uf, Utfi) 

m L J 

= — go^i + Akiqi — Bq\ + cubic order terms, (4.8) 
m 

where A and B are coefficients determined by the form of the higher-order expansion of 
e(q, u>k(q)). Note that, by Eq. (|4.7p , the leading order terms in qi in Eq. ( |4.8|) have cancelled, 
giving the following leading order behavior in ki and q\. 

de{q,u k {q)) 



-Ak 1 -2Bq 1 
oq 

(4.9) 

ql oc ±k^ 2 . 

From Eq. ( [4.2[ ) and Eq. ( [4 . 9| ) , the asymptotic form of the Born approximation plasmon 
scattering rate, as k — > k c , is 

r p (jfe) oc oc k- 112 = (k - k c y l l\ (4.10) 

Note that if up to the (n — l)th order term in q\ also cancels away in Eq. ( |4.8| ) (which 
would occur if all derivatives up to {n — l)th order of u>k c (q) and u> p (q) at q = q c are equal), 
then e(q, uj) ~ Ak + Bq n which leads to a divergence of the plasmon scattering rate which 
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goes as oc k^ 1 ^ 71 . For example, if the curvatures of the curves oo p (q) and ujk c (q) at q = q c 
were the same, then one would have a stronger divergence (k — k c )~ 2 ^ 3 in T p ^). 

This one-sided divergence in the scattering rate indicates that as the bias voltage between 
the emitter and the base, Veb, is increased so that the k of the injected electrons increases 
above k c (or alternatively, if n is decreased so that k c falls below k), the jump in the scattering 
rate should be spectacular. Thus, when V e \, is increased through this scattering threshold, 
the current passing from emitter to collector should fall dramatically. Thermal and impurity 
effects will broaden the divergence in the scattering rate, but, as can be seen in Fig. [14J this 
effect persists up to relatively high temperatures. 

A divergence in the scattering rate at the optic-phonon emission threshold has also been 
predicted in undoped one-dimension systems, due to the bandedge E~ x l 2 divergence of 
the density of final states to which the electrons are scattered. In the case we have discussed, 
the density of final states is finite (because the electrons are not scattered to the bottom of 
the band), but there is an inverse square-root (or possibly stronger) divergence in the joint 
density of states caused by the coupling of the initial and final states via plasmon emission. 
The advantage of a device based on plasmon emission over one based on optic-phonon 
emission is that the threshold energy in the former can be tuned by changing the density 
of the doping in the base region, whereas in the latter it is fixed. In GaAs, optic-phonon 
emission threshold energy is 36meV, but at present, the narrowest quantum wires have an 
energy separation between the first and second subbands on the order of about 5meV. f| 
It will take considerable technical innovation before one has good quality quantum wires in 
which the subband separation is larger than the optic phonon energy. One more problem 
with utilizing the optic-phonon principle is that the final states at the bottom of the band 
must be unoccupied, and therefore one cannot dope the system heavily, making it difficult 
to make contacts and use it as a three-terminal device. 

As the temperature is increased, the divergence becomes a finite peak. This is because 
the plasmon line broadens due to Landau damping (the collective excitations of the system 
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can have a finite life-time because they can lose energy to the single particles of the system) . 
The shift of the peak is due to an upward shift in energy of the plasmon dispersion curve 
with increasing temperature, which is a well-known phenomenon in plasma physics. [[73|] In 
one-dimension, using the Vlasov formulation [[73J] (which reduces to the classical limit of the 
RPA for equilibrium plasmas), the explicitly upward shift in the plasmon dispersion with 
increasing temperature is given, for small q, by 

m nV c (q) 

(4.11) 

<B . >= | [l + ^(W n+0(n 

where < v 2 > denotes the average v 2 over the distribution of the electron gas in the base. 
The sharp drop in the mean free path persists to relatively high temperatures (here on the 
order of tens of degrees for the parameters chosen), and therefore should be experimentally 
observable. We believe that this sharp drop in the inelastic mean free path should produce 
a large NDR in quantum wires as the injected electrons pass through the threshold energy. 

We conclude by mentioning that in higher dimensional electron systems there is a plas- 
mon threshold as well at the onset of plasmon emission. The effect in higher dimensions, 
however, is not dramatic because the ideal mean free path does not change from being in- 
finite below the threshold to zero above (as it does in the one-dimensional system) since 
single particle scattering contributes in higher dimensions, in contrast to one dimension. 
Thus, our proposed NDR in quantum wires is a specific one-dimensional many-body prop- 
erty. We also mention that a doping density of approximately 10 6 cm -1 implies that a base 
region of, say, 3000 A has only 30 electrons, which may cast doubt on the validity of our 
theory. However, various simulations of very small number of particles (on the order of 10) 
in interacting systems (for example, in the fractional quantum Hall regime or the Hubbard 
model) are sufficient to show collective effects, and therefore it is entirely plausible that the 
predicted effects will be seen even with a relatively small number of electrons in the base. 
Furthermore, with further improvements in technology, the width of the wires will decrease, 
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leading to an increase in the allowable doping density and the number of particles in the 
base region. 



V. PLASMONS IN QUANTUM WIRES 

It is apparent from the foregoing discussions that the behavior of the plasmons, especially 
in the presence of scattering and at finite temperatures, play an important role in the many- 
body physics of electrons in quantum wires. Therefore, it seems appropriate for us to 
conclude this paper with a brief discussion within the random phase approximation of the 
effects of impurities and finite temperature on the plasmon dispersion in quantum wires. 

The RPA plasmon dispersion is given by |T8| , |37| the vanishing of the dielectric function 
e(q,u) = 1 — V c (q)H Q (q, uj). The full dispersion for a pure one-dimensional system at T = 
within the RPA is [|2T],[29| 

where A(q) = exp[qn /mV c (q)]. An expansion |29|] to second order in q/kp gives, 



2 



1/2 



Wp(<7) = l?l v 2 F + -v F V c (q) + 0(q 6 ), (5.2) 

L 7]- J 

which agrees exactly with the dispersion of the elementary excitations in the Tomonaga- 
Luttinger model, a point which we return to later. For any reasonable confinement model 
V c (q — > 0) ~ |ln(ga)| and, therefore, uJ p (q) ~ \q\ | ln(ga) | 1 / 2 in one dimension. The fact 
that these one-dimensional plasmons are ungapped at uj —>■ is simply a consequence that 
the Coulomb interaction does not provide long-ranged restoring forces for charge density 
perturbations in one-dimension. 

In this section, we calculate the plasmon dispersion for a one-dimensional quantum wire 
for an impure system and at finite temperatures, by numerically finding the zeros of the 
dielectric function on the complex frequency plane, i.e., the complex frequencies u) p (q) such 
that e(q,uj p (q)) = 0. The Re^pt^)] gives the frequency of plasma oscillation, while the 
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Im[a?p(g)] gives the damping rate of the mode. As in the previous sections, we use the V(q) 
for electrons confined in an infinite square well of width a, and the RPA form of the dielectric 
function, modified by the Mermin formula to include effects of impurities. Since the modes 
oscillate as e~ iu!p ^' t , the condition of stability of the system implies that Im[u;p(g)] < 0, 
and thus we search for the zeros of e(q, z) on the lower-half complex frequency plane. Note 
that for Im[z] < 0, e(q, z) is given by the analytic continuation from the upper-half to the 
lower-half complex frequency plane. |75| Thus, the polarizability in the RPA goes as 



m 

n (g, z) 



°° dk /eq(fc + q/2) - /eq(fc ~ g/2) 

-00 k — z/q 



(h 2 )qn 

+2m(f eq (z/q + q/2) - f cq (z/q - q/2))e(-lm[z}) 



(5.3) 



where the second term in Eq. ( |5.3| ) is required for the analytic continuation for Uo(q,z) on 
the lower-half complex plane. 

As mentioned in section O, the presence of impurity scattering within the wire causes 
electrons to be diffusive at long times and large distances. This behavior of the electrons 
overdamps the plasmons and causes them to disappear at small q, as in the two dimensional 
case, J76J and this disappearance of the plasmon spectral weight at small q has important 



consequences for the the characteristics of the Fermi surface of a one-dimensional electron 
gas, as we have shown in section |I[ Here, we explicitly show in Fig. |l5|(a) this disappearance 
by calculating the plasmon dispersion for a one-dimensional quantum wire in the presence 
of impurity scattering. In the inset of Fig. [15|(a), we show the experimental one-dimensional 
plasmon dispersion compared with the RPA theory. We trapped the zeros of the dielectric 
function by using the simplex method J77|] to search for the minima of \e(q, z)\ 2 . The real 
and imaginary parts of u> p (q) are plotted above and below the x-axis, respectively. One can 
clearly see the overdamping of the plasmon mode [i.e., Ke[u p (q)] = 0] as q goes to zero at 
a critical wavevector q^. With increasing impurity scattering, q^ increases, as can be seen 
from the expansion of the dielectric function. For small q and u, the polarizability within 
the Mermin formula can be expanded to yield 
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This form can also be obtained by the standard replacement of uj 2 by u(u + 27) in the 
denominator of the usual high-frequency result for the pure system. |37j With this form of 
Eq. ( |5.4|) , the plasmon dispersion is given by 



u ( q )~^l±j-£ + ?£v( q ). (5.5) 
From Eq. ( |5.4|) , the critical at which uj becomes completely imaginary is 

qd = (5 - 6) 



where K = yme /(8ne 2 ). In the overdamped region, Im[c<;p(g)] decreases with increasing 
7 because a larger 7 suppresses the electron diffusion rate, slowing the relaxation of the 
long-wavelength density perturbations. 

We next examine the effect of finite temperature on the plasmon dispersion. Finite 
temperature effects are included by calculating U (q,u) at T / 0, using Eq. (|5.3| ) above, 
where the f e q{k) is given by the Fermi functions at temperature T. The results of our 
calculations are shown in Fig. |l5|(b) . In contrast to the case where impurity scattering 
exists, at finite temperatures, the plasmon mode at low q remains undamped. This is 
because Landau damping is caused by the transfer of energy from the collective mode to 
particles traveling at the same velocity as the (phase) velocity of the the collective mode, 
and, therefore, significant Landau damping at a wavevector q can occur only if f eq (u p (q) / q) is 
non- negligible; i. e., w p (q) / q must be less than or on the order of the thermal or Fermi velocity 
(whichever is larger). At small q, the plasmon phase velocity uj p (q)/q ~ lln(ga)! 1 / 2 — ■> 00, 
and hence Landau damping is negligible in this limit. As q increases, the phase velocity of 
the plasmon decreases and hence Landau damping increases. 

Before we conclude this section, we argue that the plasmons that we are have described 
above (which are associated with the vanishing of the dielectric function) are precisely equiv- 
alent to the Tomonaga-Luttinger bosons which are the exact elementary excitations of the 
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Tomonaga-Luttinger model. In fact, as mentioned earlier, the dispersion of the Tomonaga- 
Luttinger bosons to order q 3 is exactly the same as the plasmons calculated within the RPA, 
p9| a fact which may be attributed to the fact that only RPA-type Feynman diagrams 
(i.e. the bubble or ring diagrams) do not cancel out in the Tomonaga-Luttinger model. As 
further evidence of the equivalence of the Tomonaga-Luttinger bosons and what we call plas- 
mons, we note the similarity of the underlying physical notions of both concepts. Plasmons 
are simply long-lived charge density excitations of the one-dimensional electron gas. Since 
Tomonaga-Luttinger bosons are created theoretically by acting on the ground-state wave- 
function by charge density operators, they are also elementary excitations caused by charge 
density perturbations. Therefore, both "Tomonaga-Luttinger bosons" and "plasmons" are 
simply different names ascribed to the same physical phenomenon. This may explain why 
the experimental results || for one-dimensional plasmon dispersion (Fig. |15|(a) inset) agree 
so well with the RPA theory, which generally is regarded to be a poor approximation in 
lower dimensions. 

A similar correspondence between the spin excitations of the Tomonaga-Luttinger model 
and that of the single-particle spin density continuum can also be made. In the Tomonaga- 
Luttinger model, with no spin- dependent interactions, the spin-density excitations have a 
linear dispersion u = vfQ- The spin-density excitations (described by the imaginary part of 
the spin-spin correlation function, which, within the RPA, is exactly equivalent to the irre- 
ducible charge-density polarizability U (q,uj)) also have a dispersion of u = vpq in the limit 
that q — > 0. Unlike the single-particle charge excitations, which have a vanishing spectral 
weight as q — > for the RPA in one-dimension (all the spectral weight is pushed up to the 
plasmon [Tomonaga boson] excitations), the spin-density excitations have a nonvanishing 
spectral weight in the small q limit. This is because the spectral weight of the charge exci- 
tations is given by the screened (reducible) polarizability, Yl (q,uj)/e(q,uj), which is screened 
out at small q by the dielectric function, whereas the spectral weight of the spin excitations 
is given by the unscreened (irreducible) polarizability U (q,uj), which is not suppressed by 
the screening at small q. 
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VI. SUMMARY 



In this paper, we have studied the many-body properties of semiconductor quantum wires 
by utilizing a model that is more realistic than the exactly solvable Tomonaga-Luttinger 
model, but which we have solved approximately using the GW approximation, a renormal- 
ized perturbative scheme that has been employed with success in higher dimensions. Our 
approximate calculation provides us an opportunity to isolate the physical processes which 
are responsible for the unique properties of quasi-one-dimensional interacting systems. In 
particular, we find that the well-known result that inter-particle interactions cause the Fermi 
surface to disappear in one dimension is due to the virtual excitation of the plasmons at the 
Fermi surface. When impurity scattering is included, the Fermi surface reappears because 
the low energy virtual plasmon excitations which are responsible for the disappearance of 
the Fermi surface are suppressed by impurity scattering. This is consistent with experimen- 
tal findings which seem to indicate the existence of a sharp Fermi surface in semiconductor 
quantum wires through large Fermi-edge singularities. || Furthermore, theoretical studies 
PI indicate that for a pure Luttinger liquid which is incompressible (i.e., where the velocity 
of the plasmons goes to infinity as q goes to zero, as is the case studied here) there is no 
enhancement in the absorption cross section at the Fermi surface. [[78] 

While both our model (for zero impurity scattering) and the Tomonaga-Luttinger model 
exhibit non-Fermi liquid behavior with the interacting system having no real Fermi surface, 
there are important differences between the two models. In our model, Im[E(&F, a;)] ~ 
\u\ | lia | Co' 1 1 1 / 2 for small cu whereas in the Tomonaga-Luttinger model, Im[S(/ci?, to) ~ M", 
where the parameter a (< 1) depends on the short-range interaction constant of the model. 
The momentum distribution function n& in our model has a logarithmic inflection point 
at kF whereas in the Tomonaga-Luttinger model the inflection point in is a power law. 
We emphasize that our model is an approximate (leading-order perturbative) solution of a 
realistic model of one-dimensional electrons interacting via the long-range Coulomb inter- 
action whereas the Tomonaga-Luttinger model is an exact solution of an artificial model of 
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electrons with linear energy dispersion and infinite band width interacting via a short-range 
interaction (in particular, an interaction which is finite at q — 0). An important result of our 
calculation is that in the presence of impurity scattering the two models produce essentially 
numerically indistinguishable results. 

We have also calculated the finite-temperature self-energy of a semiconductor quantum 
wire. To perform this calculation in a numerically efficient manner, we have developed a 
formalism for directly analytically continuing the thermal Matsubara self-energy to real fre- 
quencies. Applying this method to the GW approximation, we obtain the finite-temperature 
generalization of the "line and pole" decomposition, which in which the real-frequency in- 
tegration in the spectral representation is converted into a more computationally efficient 
imaginary-frequency sum. 

In the absence of impurity scattering and at T = 0, the leading order inelastic electron- 
electron scattering rate is zero for electrons with energies below the threshold energy for 
emission of plasmons but diverges at the plasmon threshold itself. The divergence as one 
approaches the threshold from above is generally inverse square-root, and is due to the fact 
that the joint density of states of the electron states connected by the plasmon emission 
has the inverse square-root divergence (there are very special circumstances when the joint 
density of states is larger, and the divergence is in fact stronger). This sudden onset of a very 
large T t (k), which survives both the effects of impurity scattering and finite temperature, 
can, in principle, be used in a a one-dimensional THETA-type device with the possibility 
of a large negative differential resistance. This characteristic has potential applications in 
switching devices or oscillators. 

We have calculated the plasmon dispersion, including the imaginary part which describes 
the damping of the plasmons, for impure wires and at finite temperatures. Impurity scat- 
tering causes the plasmons to overdamp and disappear at small q, but finite temperature 
does not because the phase velocity of the plasmons diverges as q — > 0. 

We also calculated the bandgap renormalization at zero and finite temperatures. It is 
found to be on the order of 10 - 20 meV for typical experimental parameters, which are 
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consistent with the currently available || experimental results (which have, albeit, multiple- 
subband occupancy). We find that the band-gap renormalization can change significantly 
between T = K and 300 K for wires with very low electron densities. 
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APPENDIX A: FORMALISM FOR SELF-ENERGY Ey (*,«/„) 
FOR THE CASE OF MANY SUBBANDS 

In this appendix, we give the formalism for calculating the self-energy T,ij(k,iv n ) for 
the case where several subbands are relevant. |f{J In this case, the interaction matrix el- 
ements and the Green's functions are labeled by energy subband indices corresponding to 
the incoming and outgoing electron lines, in addition to having the momentum label for the 
x-direction. (We employ the thermal Green function formalism in this appendix.) 

When several subbands are relevant, the Coulomb matrix element has indices relating to 
the transitions from different subbands. We define Viji m (q) as the scattering matrix element 
for i — > j and I — > m, given by 

/oo /*oo 
dy d y 'v(q,y-y')UyWyMy'W m (y% (Ai) 
-oo J — oo 

where v(q,y — y') is given by Eq. (|2.3|) . 

The Dyson equation is written SIS db matrix equation 

G{k, ivn) = G°{k, %v n ) + G°(k, %v n ) E(fc, %v n ) G(k, iv n ) (A2) 

where the underline indicates the fact the the Green functions and self-energies are matrices 
and the products are matrix multiplications. The is the bare noninteracting Green 
function 
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where Ej is the subband energy. Solving Dyson's equation gives 

-l 



G 



(G ) - X . (A4) 



The quasiparticle energies are then given by the poles of G; i.e., at the frequencies such that 
det^G )- 1 -^] =0. 

The self-energy in the GW approximation [24] which is attached to external electron 
lines from subbands i and j is given by 

/dc[ 
—Wuij(q, iv m ) G° a {k + q, %v n + iv m ) (A5) 

where Wiji m =< il\W\jm > is the screened interaction between scattering states i — > j, and 
/ — > m, and the frequency summation is over the Bose frequencies iv m = ilismT . 

To calculate the components of the self-energy one must know the components of the 
screened Coulomb interaction W_. Dyson's equation for the screened interaction, for the 
random phase approximation (shown in Fig. |T|(b) ) , which translates to 

W = V + VVW, (A6) 

where V_ is the bare polarizability (given by the bubble in the random phase approximation). 
[Note that in Eq. ( |A6| ) , the matrix multiplication implies summation over two of the four 
indices, i.e., in Wiji m the "row" index is i and j, and the "column" index is m andn, and 
therefore the multiplication implied by, say, X_ = W V_, is Xij mn = J^m, W% ^vVfiumn] We 
can solve for W_ from Eq. (|A6|) , yielding 

W = e~ l V, (A7) 

where 

e = l-VV (A8) 

is the dielectric matrix. Hence, to find W_, we must invert the dielectric matrix. 

50 



The components of the bare polarizability matrix are given by 

Viji m {q,iv n ) = SuS jm n.i m (q,iu n ) (A9) 

where, in the random phase approximation, 

T-r / ■ \ fdk fi(k + q)-f m (k) 

U lm (q,iu n ) = 2 / — — — (A10) 

J 2tc £ k+q + Ei - 4 fc - E m - %v n 

where f m (k) is the equilibrium Fermi distribution in subband m. Therefore, the components 
of the dielectric matrix are 

eiji m (q,iv n ) = 5 a S jm - t^m(g)nj m (g,«V n ). (All) 

Thus, all the ingredients for obtaining W_ are available, and computing its components is 
simply a matter of algebra. 

We now give of the expressions for the Wijim(q,iu n ) for cases where (i) there is one 
relevant subband, and, (ii) there are two relevant subbands in a symmetric wire. In a 
symmetric wire, the transverse wavefunctions <j> n {y) are symmetric (antisymmetric) about 
y = for even (odd) n, and since the Coulomb potential v(q,y) is symmetric about y, Eq. 
JJ) leads to the following symmetry relationships: 



Vijim = if i + j + I + m = odd (A12) 

and 

Vijim Vjilm Vij ml Vjiml ^Imij Vmlij Vlmji Vmlji (^1^) 

(i.e., the ordering of the indices in the first two or the last two, or the ordering of the pair 
of indices, does not matter). These symmetry relations simplify computations significantly. 

In the first case, when the only one relevant subband, as in the case discussed in the 
main text, then all equations are scalar equation, and the screened interaction is (as in the 
translationally invariant 3-dimensional case) 

1 - Vu u(q) rin (q,iv n ) 
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In the second case, where there are two relevant subbands for a wire that symmetric 
about the y-axis, the dielectric matrix is 



11 

22 
12 
21 



11 

l-V nil U u 
— V1122II11 



where we have used the symmetry relations of Eqs 
Taking the inverse of e gives 

11 22 12 

^ i— V22 22II22 
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— Vn22h[22 
1 — V22 22h[22 
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1 — ^12121^12 — V12I2II21 
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1 — Vl2 12^21 

e intcr 

Vl2 12lI]J 

Q liter 
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Vl2 121^21 

e intcr 

1 — Vl2 12ni2 
C inter 



where 



(A15) 



(A16) 



(A17) 



1 - V nn (q)Tlii(q,iis n )) (l - V2222(g)n 2 2(?, iv n 

~ V ll 22(?) n il(? ; ^n)n 22 (g, Wn) 
einter(g,^n) = 1 ~ Vi 2 12(g) [Jlli(q, lV n ) +II21 (?,«!/„) 

correspond to the intra-level and inter-level dielectric functions, respectively (the zeros of 
which give the collective modes of the systems. ||21|| ) The screened interaction matrix is given 
by 

11 22 



W = e~ x V 
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22 
12 
21 
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^intra e intra 
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(A18) 
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Eq. ( |A~5| ) is then used to obtain the self-energy ^(k, iv n ). 

We now briefly discuss the effect of the second subband in the case where only the first 
subband is occupied. Then H22 = 0, and 

^ /, • n f dqfVn 11 (q)G° n (k + q, iv n + iv m ) 

E u {k,zu n ) = —T / — — r 

^ J 2tt\ 1 - Vim (q)Ilu{q,iv m ) 

+ Vi2i2(q)Gl 2 (k + q,iv n + iv m ) \ 

1 - Vi2i 2 (g)(ni 2 (g,iz/ m ) + n 2 i(g, w m )))' 

The first term in the integrand of Eq. ( |A19| ) is exactly the same as the expression we used in 
the main text to calculate the self-energy, when we ignored the higher subband contribution. 
The second term is the additional contribution of the second subband. The denominator of 
the second term is the intersubband dielectric function, [21] which only becomes significant 



for energies on the order of the intersubband energy separation. Thus, for the processes in 
which the energy is much less than the subband energy separation, as we have assumed in 
this paper, the higher subband contribution is negligible and can be ignored. 

APPENDIX B: CALCULATION FOR E(k F ,u = £ kp ) 
FOR SECOND ORDER DIAGRAM 

We show that for a short-ranged interaction, the second order diagram shown in Fig. 
|l](c) evaluated at k F and u = is identically zero. 
For V c (q) = Vq expression for the self-energy is [|T8f 



f r 1 

E(k,u) = -7^X2 J d 1j d l' ^( n F(£k+ q >)[n F (tk+ q ) - n F (€k+q+ q >)\ 

+n F (£k+ q+q >)(l ~ n F (£ k+q ))) (Bl) 



The Fermi functions restrict the four integration regions to the q, q' space shown in Fig. [16 
We calculate the contribution for each region. 
For regions I, II, and III, we obtain 

I f £ dq f £ dq' _ f 1 dx f 1 dx' 



e^oJo q Je- q q' Jo x Ji-x x' 
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f 1 , ln(l — x) 7T 2 

= — dx = — 

Jo x 6 

ii = r d -i r *l ( B2) 

JO q J2k F -q q' 

1 dx r 1 dx' 7r 2 

o x' Ji-x x' 6 

2k F dq r-i dq' r~ l dy r~v dy' 



III 



oo q J-q-2k F q' J-oo y J-y-i y' 

1 dy ( y \ f 1 , i n (i — x ) 71-2 



J_oo y v l + y 7 Jo 



ooy v l + y y Jo x 6 



(B3) 



A calculation similar to that of region III also gives a contribution of —ir /Q from region IV. 
Therefore, the second order contribution to Ti(kp,u = 0), which is given by the sum of the 
contributions from I, II, III and IV, is identically zero. 

APPENDIX C: FORM OF [T,(k F ,uj)] AS uj -> IN THE 
GW AND RPA APPROXIMATION IN ONE DIMENSION 

In this appendix, we give the details of the calculation yielding the uj —>■ forms of the 
imaginary part of the retarded self-energy in the one-dimensional electron gas, within the 
GW approximation, both with and without impurities. 

First, we define two useful energy scales 

_ k 2 F ne 2 _ 2e 2 k F 

2m e (vreo) 

Ep is the Fermi energy, and E pot is the potential energy for two electrons separated by the 

average distance between electrons at density n. As in higher dimensions, we define the 



ratio of these energies as |39] 
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TT€qVf Ep 

where v F = kpjm is the Fermi velocity. 

As mentioned previously, the only contribution to the imaginary part in the GW ap- 
proximation comes from S po i e . The contribution to Im[£] comes from two sources: (i) the 
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single particle continuum where Im[II(g, u)] 7^ 0, and from (ii) the plasmon contribution 
(i.e., when the Re[e(g, u)] = and |Im[e(g, u)\ = + ). The ^-functions in Eq. (|2.11cj) dictate 
that the only contribution to Im[E] occurs when u — £k+ q and £k+q have different signs. 
For k = k F , £,k F +q = q 2 /(2m) + qv F , and hence for small u, where one can linearize the 
energy-momentum curves about q = and q = —2k F , yielding 



Im[E(A;ir, to)] 



j / v f dq 

— Im 
2tt 



V c {q) 



+ 



dq' 

— 1m 



e(q,qv F - u)- 

V c {-2k F + q') 



jv F 2tx ^t(—2k F + q', —q'v F — 10) 



(C3) 



The bold lines in Fig. [2] show the energy-momentum curves over which the there is a con- 
tribution for the self-energy. We show below that, as uj — > 0, the Im[S] is dominated by the 
plasmon contribution in the case of clean wires, and this contribution leads to a limiting 
form o£\1m\Z(k F ,u)]\ ~ \u\ \ logd^DI 1 / 2 ). 

1. Plasmon Contribution to lm£(k F ,u) in a clean system 

The plasmon contribution comes in when the real part of e goes through zero (corre- 
sponding to the energy- momentum curve cutting through the plasmon dispersion curve). 
The contribution goes as 



Im[Z(k,u;)} = -~V c ( q *) 



9e(q,Ck+g ~ w) -1 

q=q* 



dq 



(C4) 



where q* is the wavevector at the intersection of the energy curve — lu and the plasmon 
dispersion curve u p \(q). 

For small u>, the plasmon intersection q* occurs at small q, and so one can linearize the 
energy-momentum curve 



, k F q 

t,k F + q - v ~ ^ 

m 



(C5) 



which implies that, for an intersection at q*, the corresponding energy at which it occurs is 
E* = k F q*/m. 
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The dispersion relation for the plasmon, u p \(q), is given by the solution of e(q, u p i(q)) = 
0. In the random phase approximation in a quasi-one-dimensional system, the dispersion 
relation is given by p9l 



+ — ■ (C6) 

V TIVf 

We now study two different cases, the Coulomb interaction (where V c (q) — > oo as q — > 0), 
and the short-ranged interaction (where V c (q) is finite when q — > 0). In both cases, we 
obtain non-Fermi-liquid behavior from the GW approximation to the self-energy. 

Coulomb interaction - - For the case of the Coulomb interaction where V c (q) ~ 
2e 2 | In |ga||/eo, the plasmon dispersion for small q is 



l2v F V c (q) 
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^pi(g) ~ 9\/ ~ QV F Jr s \ \n\qa\\. (C7) 



The intersection of the energy-momentum curve (,k F +q — w and the plasmon curve u) p \(q), 
from Eq. ( |C5|) and (p7|), occurs at 



q*v F (yfr s \ln\q*a\\ + l). (C8) 



In the limit of small |a>| and |g*(u;)|, the 1 on the right-hand side of Eq. ( |C8| ) can be ignored, 
yielding 

. 1 



Q 



v f Jr s \ In |g*a| 



(C9) 



Uf Jr s | In \au / (y/r~ s v ¥ )\ \ - ln(^/|ln|g*a||) 

Since \q\* <C \u\ in the limit to — > 0, the second term in the denominator of the last expression 
in Eq. (|C9|) is negligible when compared to the first term and hence 

?» « - 1 (CIO) 

^/r s | ln|aa;/(^u F )|| 
Now the derivative of the dielectric function for q — > and u/q ^> v F is 
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de(q,£k+q-w) dV c (q)U(q,^ k+q 



to 



dq 



dq 

n (q, u) + V c (q) 

dq m dq 

n dV c (q) q 2 

m dq (qkr/m — u) 2 

2q 



+ 71 y (g) d (<l 2 / (l k F/m - u) 2 ) 



+ ln(M) 



m dq 
2(q 2 k F /m) 



(Cll) 



\qkpjm — oo) 2 V (qkp/m — uj) 2 (qkr/m — lu) 3 ' 

Given that q ~ oj/\ ln(|o;|)| <C uj as uj — > 0, the second term in the last expression of Eq. 
( |C11| ) dominates, and thus we have 



9e(q, £ k+q - w) 



dq 



<?=<? 



r s w F | ln|g a||— . 



(C12) 



Substituting Eqs. ( PT2|) and (|CT0|) into Eq. gives 



Im[£(&F , u;)] 



7T\/T 





In 















(C13) 



Thus we have shown that in the clean limit, within the GW and RPA approximations the 
imaginary part of the self-energy goes as 



|Im[£(/cF, oj)]\ ~ In |u;| |, 



(C14) 



which implies there is no Fermi surface, and we have shown that the physical mechanism 
for the disappearance of the Fermi surface is the virtual emission of low energy plasmons. 

Short-ranged interaction — In the case of a short-ranged interaction (i.e., & delta-function 
in real space) where V c (q) = U is a constant, we show below that |Im[£(/cp,u;)]| ~ |a>| for 
small uj. 

The plasmon dispersion for short-ranged interactions is 



^ P i(<?) = qVF\jl + 



TlVp 



(C15) 



The dielectric function at small q for z = — w is given by 



k F +q 



UJ 



mUo 



In 



nq 



1 + 



VFq 



m(uj 2 — 2u>vpq) 



1 - 



vpq" 



m(uj 2 —2u)Vpq) 



1 - 



2U v F q 2 



tt(uj 2 — 2uJv F q) 
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(C16) 



and therefore, 

de(q,^ kF+q - uj) 2U v F [ 2q 2uJv F q 2 



(C17) 



dq 7r ^{uj 2 — 2uJVpq) (uj 2 — 2ujVpq) 2 - 

Substituting Eq. ( pl7| ), into Eq. ( |(J4j) yields the plasmon contribution to the self-energy, in 



the limit of to — > 



tt(,/1 + 2U /(irv F ) - Ifuj 
Im[E pl (A*,o,)] = - W . (C18) 

Ul + 2U /(m) F ) 



Thus in the case where the interaction is short-ranged, one obtains the imaginary part 
of the self-energy being proportional to uj, giving a marginal Fermi liquid. From exact 
solutions, we know that short-ranged interactions give non-Fermi liquids. We have shown 
that the GW approximation is sufficient to describe this non- Fermi-liquid behavior for short- 
ranged interactions. (Note, however, that the functional dependence is different in the exact 
Luttinger liquid situations, where Im[£(/cp, lu)] ~ \oj\ a , with a < 1.) 

2. Single particle contribution to Im[X(&F,a;)] in a clean system 

The single-particle excitation contribution to Im[E] comes from each of the two terms 
on the right-hand-side of Eq. ( |C3| ). The first term corresponds to virtual electron-hole 
excitations around q = 0. Within the region of integration < q < uj/vp, Im[e] 7^ only in 
a region the size of ~ uj 2 . In this region, both the real and imaginary parts of e(q, qvp — uj) 
are of the order of V c (ui/vf)/uj, and therefore the integrand V c (q) Im[e _1 (g, qvp — lS)\ is on 
the order of u. Thus, the contribution of this first term is on the order of the product of the 
size of region of integration and the magnitude of the integrand, uj 2 x uj = uj 3 . 

The contribution from the second term of Eq. ( |C3|) corresponds to virtual electron-hole 
excitations around q = 2kp. Here, the region where lm[e(2kp + q', —q'vp — uo)} 7^ is of order 
uo. In this region, Im[e] ~ constant, and Re[e] ^ | ln(o;)|, leading to the integrand being on 
the order of |ln(a;)| 2 . Therefore the magnitude of the contribution from this term is the 
product of the size of the region of the integration and the magnitude of the integrand, 
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which is uj/\ ln(u;)| 2 . Note that the ln|u;| comes from the form of the RPA Uo(q,uj) and 
is independent of the form of V c (q). This term dominates the q = term, and hence the 
single-particle contribution to the self-energy is 

\I m [Z(k F ,u;)}\~-M—. (C19) 



3. Im[£(/cF, lo)] as lo — > for dirty systems 

In the case where there is impurity scattering, electrons are diffusive and, as shown in 
section [V], the plasmons no longer exist at small wavevectors. Thus, for small to, there is no 
plasmon contribution to the imaginary part of the self-energy. In this subsection, we show 
that the remaining contribution goes as Im[£(fcp, u;)] ~ lo 2 . 

We find it convenient to define the following the following dimensionless variables: 

k n k q o u 

kp kp hp 

7 = V(Q) = ^T~, (C20) 

E F 2e 2 /e 

n(Q^) = 7 H #T, m,n) E ^ u) nk ' u) 



(2m/Tck F )' "~ E pot 2e 2 k F /(ire ) ' 

With these definitions, fl(Q -»• 0, Q = 0) = -1 and e(Q, Q) = 1 - r s V c (Q)fl(Q, ft). In terms 
of these dimensionless variables, Eq. (|C3|) becomes 

dQ ~ 



rU A (It J _ 

E{K = 1, ft) « / K(Q) Me (G> 2 Q " «)] 

jo 2 



+ / Vc{-2 + Q') Im[ e - 1 (-2 + Q', -2Q' - ft)]. (C21) 

J-n/2 2 

The polarizability is given by the Mermin form. In the region Q, O C 7, the polariz- 
ability is of the form 

■W^-Q^Ti <C22) 

which gives 
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1 = 1 tfjg . (C23) 

Therefore, the virtual transitions around Q = [the first term in Eq. (|C21j )] give a contri- 
bution 

Tsl Z^ 2 Q 2 (Q-2Q)V 2 (Q) 



U ' ^ 8 io V g4 (1 + rsK( g ))2 + ((fi _2Q)^/4)2 



(C24) 



27 7o (l- y )4((i + rs y c (^(l-y)/2))0/7) 2 + y 2 ' 
As Q — > 0, the integral in the last line of Eq. ( |C24|) diverges in the region y = 0, and 

therefore only the behavior of the integrand around that region is relevant. Therefore, the 
asymptotic form of S (using the asymptotic form V C (Q) ~ | ln(Q)|) is 

r s 7 /"! hi(O) 2 ?/ 



E 1 (A' = l,fi)«^ / rfy 
o Jo 



[uj 2 r 2 ln(o;) 2 /7 2 ) + y 



2 



rtfMar^Mary (C25) 



27 V 7 

which to leading order, goes as f2 2 | In |0|| 3 . 

Now we study the contribution from the Q = 2 virtual transitions. The polarizability in 
this region is 



(C26) 



ri 7 (Q, a) = n (Q, n = o) - zC(q)q + o(n 2 ), 
c (no(Q,»7)-no(g,o))n (Q,o) 

7n (Q,i7) 

Using this expansion, we obtain 

Im[e-(g, n)] - - Im[ ^ (g ^ )] , (C27) 

L 1 ;J (r s K(Q)Re[n (g,o)]) 2 ' 1 ; 

which yields a contribution to Ti(K = 1,0), from Eq. (|C21|) , of 

= fi) ~ f-«n dQ ~ !W (C28) 

Thus, in this case where electrons diffuse due to impurity scattering, the virtual excita- 
tions close to Q = dominate, and 

\ImZ(k F ,tu)\~uj 2 \ln\uj\\ 3 . (C29) 
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4. Form of Re[£(/cF, u)] for clean systems 



The u — > limit of the real part of the self-energy at k = kp is 

r n/2 aq _, 

Re[S(fc F ,w)]^ / -^^(QjRete-^Q^Q-n)-!] 
Jo 2 

+ / -T- ^(-2 + Re[e- 1 (-2 + Q', -2Q' - Q) - 1]. (C30) 

7-n/2 2 

The real part of e _1 (Q,fi) in the integration region Q = 0,Q/2 is well-behaved (aside from 
an inconsequential principal part divergence at the plasmon frequency). The e _1 — 1 term 
basically is a constant, and the magnitude of the integrand is determined by the V C (Q) term. 
Since the size of region of integration goes as Q and the integrand goes approximately as 
V C (Q) ~ | In the contribution from the first integral on the right-hand side of Eq. ( |C30|) 
is on the order of \Q In \ Q\ \. In the second integral, the integrand does not diverge, and hence 
its contribution is on the order of Q and can be neglected in comparison to the first. Hence 

Re{E{k F ,u)} ~ |wln|a;||. (C31) 

The slope of Re[E(/c^, a;)] at u = is infinite (in fact, this is necessary for A(kF,u) to be 
integrable), implying from Eq. (|2.20| ) that Zp = 0, which is consistent with the fact that 



there is no Fermi surface. 

APPENDIX D: Im[E(k F , u)] AS u -> 0, IN THE 
GW AND RPA APPROXIMATION IN TWO DIMENSIONS 

In this appendix, we calculate, for the sake of comparison, the imaginary part of the 
T:(k F , uj), due to plasmon and single particle emission in two dimensions for the case of 



Coulomb interaction, which is 26 



2-7TP 2 

VM) = — ■ (Dl) 

In the GW approximation for a spherically symmetric parabolic band, the imaginary 
part of the self-energy can be written as (for the case of u > 0) 
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2m 



lm[E(k,u)\ 
where e(q, uS) = 1 — V c (q)U(q, uj) and 



(2tt)- jo 



JO JO 



g 2 (0,£) = k 2 + kp + 2m£ - 2k J k% + 2m 75 cos 9. 



(D3) 



The exact form of U(q,uj) in RPA can be found in references |81[] and |26[. We define the 
reduced variables 



k q 

kp kp 
2mi i 



y 



kp 



V C (Q) 



Vc(q) 



Ep' 
1 



E F ' 



7 

2e 2 



(2ney(q F e ) Q'' 
u;) 



7 : 



(D4) 



n(Q,u;) 



n(g,u;) 



e 2 kp/('Keoy * m/7r 

With these definitions, S(Q -»• 0, = 0) = -1 and e(Q, ft) = 1 - R S V C (Q)U(Q, Q). For the 
case k = kp, Q 2 = 2 + y — 2y/l + y cos 6 1 , which yields 

QdQ 



dd 



sm9(Q,y) y/T+j}' 



where, 



sin 9 (Q,y) 



4 Q2 _ y 2 + 2 yQ2 _ Q4 



4(1 + y) 

Therefore, using Eqs. ( |D4j ) , ( |D5| ) and (D), Eq. ( |D2| ) can be rewritten as 

lm[t{kp,u)} = Fdy f^ V) dQlm[e-\Q,y-n)]-== * 

Jo JQ-(y) V4Q^ - JT + 2yQ 2 - Q 4 

where, 



(D5) 



(D6) 



(D7) 



Q_(y) = [2 + y - 2^/l + y] 1 = | + 0(y 2 ) 



(D8) 



Q + (y) = [2 + y + 2^/l + y] 5 = 2 + | + 0(y 2 ). 



Within the GW approximation for two dimensions, Eq. ( |D7| ) is an exact result. Now, 
we calculate the plasmon and single-particle excitation contributions to Y>(kp,uj) in the 
uj/Ep < 1 limit. 
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1. Plasmon contribution in two dimensions 



The plasmon contribution to the imaginary part of E(fc, u) occurs when Re[e(Q, y — Q)] = 
in Eq. ( [D7| ) . Since we are interested in small frequencies, the dielectric function in two- 
dimensions is given in the small frequency limit by |J7J 

««■»-«> c») 

and therefore, the wavevector at which Re[e] = is 

QM = (dio) 

Hence, 

nS(Q-QJy) 
lm[e-\Q,y-Q)\ K " ] 



\de(Q,y-n)/dQ\ 



^(Q-Q P (y)y\—L- (Dll) 



2R_ 

There is a non-zero plasmon contribution to the integral in Eq. ([D7]) only if the delta function 
peaks between Q-(y) and Q + (y); that is (since Q p (y) is small) only if Q p (y) > Q-(y). This 
occurs if the condition 

(g - y? 

R s 

is satisfied, or equivalently, if 



> y, (Di2) 



„< w = ^( 1 + __^ 1 + _) = _ + ( fl . ) . (D13) 

Thus, from Eqs. (|D7|), ( piO| ), ( |L) 1 1| ) and ( |L)13| ) , the plasmon contribution, to lowest order 
in Q 2 , is 

e 2 k F r^/Rs (y-n) 2 



Im[Z pl (k F , u)\ « — — / dy 

L€nii<, JO 



o^sJo ^(n 2 /R s ) 2 -y 



2 Jo VT 

71 U) 2 
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(D14) 



Therefore, the plasmon contribution to the imaginary part of the self-energy goes as as u 2 . 

Note that the quantity we calculate here, Im[£(fcp, u)}, is different from the inverse 
inelastic lifetime of the quasiparticle, r~ 1 (fc) = 2 ImfE^^, £k F )]- F° r Tec 1 there is a minimum 
critical \p — Pf\ before the plasmon contributes to r" 1 , ^6] whereas for Im[£(fcp, lj)], there 
is no gap. 



2. The single-particle excitation contribution 

We show in this subsection that the single-particle excitation contribution to the imagi- 
nary part of the two-dimensional self-energy goes as uj 2 \ ln(|u;|)|, and hence gives the dom- 
inant contribution to the Im[S] in a two-dimensional electron gas. This result has been 
obtained previously, 0,0] in similar contexts, although details of the actual calculation 
were never given. Below, we give a quick derivation of this result. 

For small Q, the polarizability goes as |81|J26|] 



Re{IL(Q,Q)} 



Im[n(Q,fi)] ! ^ 



(D15) 



and hence, since |Im[n(Q,fi)]| <C \Re[H(Q, in this limit, 

Me-HQ m - Im[e(g '" )] - ™ (D16) 

The single-particle contribution to the imaginary part of the self-energy, from Eq. (|D8|), is 

[mP , p( ^)] « $ jf* ( " - S) Q(Q + R%-Q>) (D17) 

where terms higher order in y and have been ignored. The lower limit for the Q-integration 
is given by a cutoff on the order of Q, which is where the expansion given in Eq. ( pi6| ) 
starts to fail (the remaining part of the integral is from ~ Q to Q- is well behaved and gives 
negligible contribution). The leading order contribution to the Eq. ( [01 7| ) is 
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lm[Z(k F , u)] w ^ /" dx (Q - x) ln(fi) 
7ren Jo 



7re j o AR S 

\n\w\ (D18) 



8vr£ F 

Thus, the single-particle contribution to the imaginary part of the self-energy is on the order 
of u 2 \ In and in the small \u\ limit, is the dominant term for Im[£(/cp, us)]. 

3. Im[E](/cp, u) for an impure system 

In this subsection, we show that in a dirty system, the main contribution to lm[T,(kF, u)] 
comes from the small Q and y region in Eq. ([07]), which gives a linear contribution in u. 
Defining x = (Q — Q-(y))/y, we rewrite Eq. ( |JJ7| ) in the limit of small Q as 



n roc Im 



l + 2x 



Im[Z(k F ,n)\ « / dyy dx — t v - z n '-± (D19) 



/o Jo |/(:r + !) 4a/^ + x~ 

We first calculate the inner integral of the right hand side of Eq. ( |D19| ) , 



OC 



Im 



-i(y(x + l),y-n) 



l + 2x 



J (y) = / dx 1 v , , ' ^-==. (D20) 



'o 2/(a;+!) 4\/z + x z 

We approximate the dielectric function e(Q, O) by using the diffusive form of the polarizabil- 
ity, written as an expansion in powers of Q and Q of the numerator and denominator of the 
Mermin formula |JT]. The diffusive form is strictly valid only for Q, Q << 7. However, since 
the contribution to Im[E(&p,0)] from the small Q region goes as Q, and the contribution 
from the large Q region goes as Q 2 both in the diffusive and the ballistic (i.e., 7 = 0) forms, 
the error introduced by using the diffusive form at large Q is insignificant. The diffusive 
form of the polarizability is 

Q 2 



n 7 (Q,n') 



Q 2 - 



Re W'> * - iyrfb? (D21) 



Imn 7 (Q,fi') 

and hence 



2(Q 4 + (7fi/2) 2 
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Q 



i? s Im[n] 



(Q - R s Re[n\y + {R s im[n\y 



2R S 
7Q 



(D22) 



where we have ignored terms higher order in Q,Q. The inner integral is therefore 

2x + 1 



i(y) 



dx 



2^ Jo ' (l + A y (2x + l) 2 )Vx 2 + x 



By a change of variables z 2 = x 2 + x, we obtain 



lyj ~ 7^ Jo 1 + A,(l 



(1 + Az 2 ) 



AnJA y (A y + l)7fi 



and hence 



JO 

1 // n 



(y - ny 



AnQ Jo ^(r s y 2 /^) + (y - Of 



(D23) 



(D24) 



, . (D25) 

/o 4vr + (1 - C) 2 

The integral is just a function of r s /7, and therefore this shows that |Im[E(/Cf,a;)]| ~ in 
a dirty system. 



APPENDIX E: ESTIMATION OF THE COUPLING PARAMETER g 

FOR A QUANTUM WIRE 

In this appendix, we estimate the strength of the Tomonaga-Luttinger coupling param- 
eter g within a quantum wire of width a which is screened by another quantum wire which 
is adjacent and parallel to it, a distance / away. We denote the wires by numbers 1 and 
2. We assume that the electrons do not tunnel and are screened within the random phase 
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approximation. The effective intrawire interaction (i.e., electrons within the same wire) is 
given by the the bare interaction screened by the motion of the electrons in the adjacent 
wire. The Dyson's equation for the screened interaction is 

V eS (q,u) = V c>11 (q) + V cA2 (q)U 2 (q,u;)V c>21 (q) 

+V c , 12 (q)U 2 (q,u)V Ci22 (q)Il 2 (q,u)V Ci21 + - 

oo 

= V c>n {q) + V c>12 U 2 {q,u){q)V c , 12 {q) £ IL 2 {q, u)V c>22 {q) (El) 

n=0 

_ T/ / \ , K,i2(g)n 2 (g,a;)K,2i(g) 
1 - V c>22 (q)U 2 (q,u) 

where V c ij is the bare Coulomb interaction between electrons in wires i and j, and ITj is the 
polarizability of wire i. 

Now assume both wires are identical, so that the intrawire interaction Vn = V 22 = V a , 
the interwire interaction Vi 2 = V 2 \ = V e and Iln = II22 = n. Since we are interested in the 
physics close to the Fermi surface, we take the q — > asymptotic forms of these potentials. 
The intrawire potential V a is given by Eq. ( |2.5| ), and the interwire potential V e is given, in 
the limit of small q, by 



2 e 2 fa/2 rl+a/2 

— / dy 

Co J-a/2 Jl-a/2 

K (\ql\) « V a (q) - 2eV ln(//a). (E2) 



V e (q) = — f 2 dy f +a/2 dy'K (\q(y-y')\) cos 2 (^) cos 2 

6n J-a/2 Jl-a/2 K a ' K d 



This gives a long wavelength effective intrawire interaction of 

V a (q)-2e% 1 \n(l/a)(2V a (q)-^% 1 \n(l/a))u(q,uj) 

V eS (q,u) = — ^— — — . E3 

1 - U(q,u)V a (q) 



In the limit q — > 0, V a (q) — > 00 and U(q,u = 0) = —(2/tyvf), which yields, from Eq. (|E3|) , 

V cS {q - 0, 0) = ^(l + ln(//a)) . (E4) 

For //a = 10, and for the parameters of Fig. |], 8e 2 jie^Vpii) = 2.8, we obtain g = 
2V e s/(7rv F ) = (l + 2.81n(10)) w 7.5. and therefore, using Eq. ( gT34|) , a w 0.15. 
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APPENDIX F: Im[E(k F ,u)] AS u -> 0, 
IN THE GW^r AND RPA APPROXIMATION 



In the GWT approximation, ]8^,iT| the Matsubara expression for the self-energy is 



p J 2tt e(q, %v m ) %v n + %v m - £ k+q 

where 

e(q, iu n ) = 1 - V c (q)Il(q, Lu)T(q, iv n ), 



(F2) 



r(q,iu n ) = l/[l + G(q)V c (q)U(q,u)). 



T gives contributions from the vertex corrections. In the GW and RPA approximations 
we have used, T = 1. We now show that when vertex corrections in the form given in Eq. 
( |F2|) are included, the main conclusions of our calculations regarding the existence of Fermi 
surface do not change. 

Substituting Eq. ( |F2"|) in Eq. (|F1|), one can rewrite the expression for the self-energy as 



" P h J 27r 1 - Vc(qMq, lj) [1 - G{q)\ iv n + iv m - £ k+q ' 1 
Physically, this equation says that the main effect that the vertex corrections have is that the 
local-field factor G(q) (which accounts for the short-range correlation of electrons beyond 
RPA) is to renormalize the dielectric function, changing the dispersion at large q. 

One can make the equivalent of the Hubbard approximation |83],f34| for G(q) in one 
dimension to obtain |85[ 



i VcUq 2 + kl) 

In the limit where q — > 0, G(q) ~ | In |ga| | _1 — > 0, and therefore, the local-field correction 
term in Eq. (|F3|) is negligible for small q. However, the behavior of Im[£(/cp, u)} for small 
uj only depends on the integrand of Eq. flFll ) at small q, which the vertex corrections leave 
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essentially unchanged. Thus the inclusion of vertex corrections of the form Eq. (|F2|) does 
not change at small u, and hence does not change the conclusions regarding the 

Fermi surface. 



APPENDIX G: EXPRESSION FOR THE FINITE-TEMPERATURE 
SECOND-ORDER SELF-ENERGY 

Here, we give the expression obtained for the second order (in the screened interaction 
w) self-energy (Fig. 0(c)) using our contour deformation method. We obtained 

=<*•*> = / (5^/(1^™ 'M (G1) 

where 

Hk / z \ T 2y^y w(q',iu n >)w(q,iu n ) 

k ' q ' q £^ t? n , fan + Z- £ k+q ) (iu n + %UJ n > + Z- £k+q+q') fan' + Z- £k+q') 

| T ^ w{q,iuj n )w{q',^ + ^ - z)R z (^ +q >) 

iuj n fan + ^k+q' - Ck+q+q'X^n + z ~ £k+q) 

+ r ^ w(q',iuj n )w(q,£ k+(l - z)R z (£ k+(l ) 

iw n fan + £k+q - 6c+q+q')(^n + Z~ £k+q') 

_ r ^ w(g',^ n )w(g,(k+ q +q' - z - ^rJ-Rz^k+g+gO 

^ ' ' ' + Z - £k+q.')fa + £k+q - £k+q+q'^ 



[lLU n 



+ w(g,£k+ q - z)w(g',£ k +q' - z)R z (^ +tl )R z (^ +fl >) 

£k+q + £k+q' _ £k+q+q' — z 
_ w(g',£k+q' ~ z)H/(z,£ k +q+q' ~ Ck+qO^ Z (^k+qO^ F (^k+q+qO 
£k+q + £k+q' _ £k+q+q' — Z 
_ f°° B ^ Ck+q+q' -Z- ^)# 2 (£ k +q+q' ~ u) [^(Ck+q+qQ + n B {u)] 

J-oo 2n (U + £k+q' — £k+q+q') (id + Z — £k+q) 

where R Z (E) = n B (E - z) + n F (E). 
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TABLES 

TABLE I. The small forms of the real and imaginary parts of S in one and two dimensions, 
in the GW and RPA approximations, for clean and dirty systems 





d= 1 


d = 2 


Im[S] , clean, (plasmon contribution) 


uj\J ln( uj ) 


UJ 2 


Im[S] , clean, (single-particle contrib.) 


^/[m(M)] 2 


" 2 \H\uj\)\ 


Re[£], clean 


uj \n{yj) 


UJ 


Im[E] , dirty 


^ 2 |[ln(k|)]| 3 


\uj\ 


Re[E], dirty 


UJ 


> UJ 



TABLE II. Many-body characteristics of one-dimensional quantum wire with and without scat- 
tering, in the random phase approximation. 





7 = 


7/0 


Fermi liquid? 


No (Z F = 0) 


Yes (Z F + 0) 


Re[£(fcp,u;)] as uj —>■ 


~ a;log(|u;|) 


~ — UJ 


|Im[£(fci?,u;)]| as uj — > 


~ |^| | logdu;!)! 1 / 2 


~ ^ 2 |log(|o;|)| 3 


A(kp,oj) as uj — ► 


~ | lo g (|^|)|) 1 


~ |log(|^|)| 3 - 


Inelastic scattering rate L 


{k-kc)- 1 / 2 ask^k+ 


Finite for all k 


below plasmon emission threshold 


Nonzero except at k F 
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FIGURES 

FIG. 1. (a) Electron self-energy in leading order in the effective dynamical screened Coulomb 
interaction. The straight line denotes the bare Green's function and the thick wavy line denotes 
the dynamically screened interaction, (b) Dyson's equation for the screened interaction, within 
the random phase approximation (the thin wavy line is the bare Coulomb interaction). Impurity 
effects are introduced in the electron-hole polarization bubble through the Mermin formula, (c) 
Second order diagram generally ignored in this paper. In appendix B, we show that this diagram, 
for k = kp and u = —£k F with a screened interaction which is a delta- function in real-space (i.e. a 
constant in momentum space), is identically zero. 

FIG. 2. Regions of q— uj' space with non-zero lm[e~ 1 (q, uj')]. The grey shaded regions correspond 
to the single-particle excitations, Im[e] 7^ 0, and the thin line corresponds to the plasmon dispersion 
Re[e] = 0. To evaluate Im[E(A;,u;)] within the GW approximation, one integrates the function 
V(q)Im.[e (q,uj')] over the two bold line segments, which belong to the curve uj'(q) = uj — Ck+q) 
[see Eq. Q2.11c| )]. The figure shows that both the single-particle excitation region and the plasmon 
line contribute to Im[E(&F, uj)]. 

FIG. 3. (a),(b) Real (bold lines) and imaginary (thin lines) parts of the self-energy T,(k,uj) and 
(c), (d) spectral functions A(k,uj) as functions of the energy uj, for 7 = (solid lines) and 7 = Ep 
(dashed lines). Figs, (a) and (c) are for k = 0, and (b) and (d) are for k = kp. The parameters 
used are kpa = 0.9 and r s = (4m e e 2 /irh 2 kpeo) = 1.4, |3J| (corresponding to a = 100 A, a density of 
0.56 x 10 6 cnr 1 and a Fermi energy of Ep ~ 4.4 MeV in GaAs). The vertical arrows in (c) denote 
5-functions in the spectral function at uj = —4.9, 0.9 with weights (2ir) x 0.32 and (2tt) x 0.33, 
respectively. Note that in (d), A(kp,u}) has a 5-function of weight (2ir) x 0.3 at uj = for 7 = Ep, 
but not for 7 = 0. The straight lines are given by u — — fj,, and their intersections with Re[E] 
indicate the solutions to Dyson's equation and correspond to a quasiparticle peak. 
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FIG. 4. (a) Momentum distribution function n/-, around k/kF = 1, (b) Fermi surface renormal- 
ization factor Zp, and (c) density of states p(ui) for small u, of a quasi-one dimensional electron 
gas for kpa = 0.9 and r s = 1.4 (as in Fig. |3|), for various impurity scattering rates 7. The bold 
lines refer to k > kp, and the thin lines to k < kp. For 7 = 0, is continuous at k = kp, implying 
that the system is non- Fermi liquid, but for a discontinuity occurs at k = kp, signalling the 

presence of a Fermi surface. The magnitude of the discontinuity in nk F is given by Zp. 

FIG. 5. (a) Inelastic scattering rates T(k), and (b) the corresponding mean free path 
l(k) = T(k)k/m, as a function of k, for various 7's (electron-impurity scattering rates), for 
kpa = 0.9 and r s = 1.4. Within RPA, for 7 = 0, the T(k) is identically zero below k = k c because 
energy and momentum conservation prohibits single particle excitations and plasmon emissions. 
Above k c , the scattering rate is caused by plasmon emissions. For 7 / 0, the plasmon line broadens 
and momentum conservation is relaxed, resulting in a nonzero T for k < k c . 

FIG. 6. Total bandgap renormalization (Re[S e + S^] at k = 0, u = £k=o) as a function of 
electron density in the quantum wire for various wire widths with parameters corresponding to 
GaAs. 

FIG. 7. Electron self-energy for impurity scattering, to lowest nontrivial order in the impurity 
potential, used to calculate the momentum distribution function of a disordered one-dimensional 
electron gas (see Fig. |8|). The dashed lines are the screened impurity scattering potential, which 
we approximate by 5- functions in real space. 
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FIG. 8. (a) The momentum distribution function, for 7i m p/-E'F = 0.02,0.2 and 2, and (b) 
spectral function (bold solid line), and real (thin solid line) and imaginary (dashed line) parts 
of the self-energy and for k = /cp, 7i mp = 2Ep, of a disordered non- interacting one-dimensional 
electron gas at T = 0, calculated with the self-energy shown in Fig. [/]. 7i mp = 2mN i\J§ j hp is 
the Born approximation impurity scattering rate for the electron of energy Ep (the Fermi energy 
corresponding to the pure system of the same electron density). The inset in Fig. (a) shows 
the chemical potential [i for fixed density, as a function of the disorder. In Fig. (b), the arrows 
correspond to delta functions. The dashed arrow marks the position of the delta function for the 
pure system, and has weight 2ir. The solid arrow has weight (2ir) x 0.39. The straight dotted line 
is given by u — — fi, and its intersections with Rep] correspond to a quasiparticle peaks. 

FIG. 9. The distribution function of a Luttinger liquid with Fermi surface exponent a = 0.15, 
which corresponds to a dimensionless coupling constant g ~ 7.5, for various scattering rates. The 
Luttinger model linearizes the dispersion around the Fermi surface, and therefore the curves strictly 
should only be valid around k = kp. 

FIG. 10. The contour of integration C for Eq. ( |3.1C| ). The hatched real axis indicates a branch 
cut due to w(c[, ui) in the integrand of Eq. (|3.10[) . The crosses mark the poles due to the integrand; 
the ones on the imaginary axis are due to ns(w), and the isolated pole is due to the denominator. 
The residues of the poles on the imaginary axis give h^^z), while the residue of the isolated pole 
gives h kt q(z). 
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FIG. 11. (a), (b) Real and (c), (d) imaginary parts of the self-energies and (e), (f) spectral func- 
tions, as a function of frequency, for electrons in a quantum wire, for k = 0, k F , at finite tempera- 
tures. As in previous figures, /c^a = 0.9 and r s = 1.4, which correspond to Ep k, 4.4 meV k, 50 K 
in GaAs. The discontinuities is Im[S(/c, u>)] at T = 0, which arise from virtual plasmon emis- 
sion thresholds, broaden with increasing temperature because the plasmon peaks broaden due to 
Landau damping. The logarithmic divergence in the T ^ ImS(/c,cj) at to = E^ — fio(T), due 
to the divergence of the one-dimensional Born-approximation scattering rate, shifts slightly with 
respect to temperature because of the temperature dependence of /io, the non- interacting chemical 
potential. 

FIG. 12. Temperature dependence of the band-gap renormalization due to conduction electrons 
for a wire width of 100 A in GaAs, for electron densities of 10 4 cm -1 (solid lines), 10 5 cm" 1 (dotted 
lines) and 10 6 cm -1 (dashed lines). The thin lines are for the electrons (Re[£ olcctron (/c = 0,^ = o)]) 
the light bold lines are for the holes (Re[£ hole (/e = 0, £fc=o)] and the heavy bold lines are for the sum 
of the two. The densities n = 10 4 cm" 1 , 10 5 cm" 1 and 10 6 cm" 1 correspond to Fermi temperatures 
of E F = 1.6 x 10" 2 K, 1.6 K and 160 K, respectively. 



82 



FIG. 13. (a) A schematic of the band diagram of a one-dimensional tunneling hot electron 
transistor, where electrons are injected from the emitter into the base (which contains a Fermi 
sea of electrons) and (b) energy- vs. momentum-loss diagram for the injected electron. In the 
transistor, the electrons that are scattered below the base-emitter barrier do not reach the collector, 
and therefore the fraction of the injected electrons reaching the collector depends on the electron 
scattering rate in the base region. In Fig. (b), the solid (dashed) line indicate the energy- vs. 
momentum-loss for electrons injected into the base region below (above) the plasmon emission 
threshold (i.e., the solid line is for k < k c and the dashed line is for k > k c ). The intersections 
of the energy- vs. momentum-loss curve and the plasmon dispersion curve (bold line) indicates 
the wavevectors at which plasmons are emitted; if there is no intersection (as with the solid line), 
plasmon emissions are not allowed. As the energies of the injected electrons is raised above the 
plasmon emission threshold, the scattering rate increases dramatically (see Fig. |l4|), drastically 
reducing the fraction of injected electrons that reach the collector. 

FIG. 14. (a) Momentum scattering rate T t ^ and (b) the corresponding mean free path, 
h = v k/^t,k- of an electron in a doped one-dimensional quantum wire, as a function electron 
momentum, for various temperatures. As in previous figures, fcpa = 0.9 and r s = 1.4. 

FIG. 15. Calculated one-dimensional plasmon dispersion in a GaAs quantum wire within RPA 
for (a) T = and various 7's, (b) 7 = and various temperatures, calculated by finding the 
zeros of e(q,u) on the complex frequency plane. The curves with uj > give the real part and 
those with oj < give the imaginary part. The curves include full V c (q) for infinite square well 
confinement. As in previous figures, the system parameters are kpa = 0.9 and r s = 1.4, which, for 
GaAs, correspond to a Fermi wavevector of kp = 0.88 x 10 6 cm -1 . In (b), we show both the long 
wavelength result, uj 2 = nq 2 V c (q)/m, using V c (q) = 2e 2 | ln(go)|/« (dotted line) and the full V c (q) 
(long dashed lines) . Experimental results of reference || compared with RPA theory (of reference 
PHI ) are shown as an inset in Fig. (a). 
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FIG. 16. Regions of integration for the calculation of the second-order self-energy, for k 
and lo = 0. 
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